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Abstract. We present the results of high angular resolution millimeter observations of gas and dust toward 
G31.41+0.31 and G24.78+0.08, two high- mass star forming regions where four rotating massive toroids have 
been previously detected by Beltran et al. (121)1)41 . The CH 3 CN (12-11) emission of the toroids in G31.41+0.31 
and core Al in G24.78+0.08 has been modeled assuming that it arises from a disk-like structure seen edge-on, 
' with a radial velocity field. For G31.41+0.31 the model properly fits the data for a velocity v rot ~ 1.7 km s _1 

at the outer radius R ou t — 13400 AU and an inner radius 7?i nn — 1340 AU, while for core Al in G24.78+0.08 
the best fit is obtained for v TO t — 2.0 km s" 1 at R ou t — 7700 AU and i?i nn — 2300 AU. Unlike the rotating disks 
detected around less luminous stars, these toroids are not undergoing Keplerian rotation. From the modeling 
itself, however, it is not possible to distinguish between constant rotation or constant angular velocity, since both 
velocity fields suitably fit the data. The best fit models have been computed adopting a temperature gradient of 
, the type T oc i?~ 3 ^ 4 , with a temperature at the outer radius T out ~ 100 K for both cores. The Md yn needed for 

equilibrium derived from the models is much smaller than the mass of the cores, suggesting that such toroids are 
unstable and undergoing gravitational collapse. The collapse is also supported by the CH3 13 CN or CH3CN line 
width measured in the cores, which increases toward the center of the toroids. The estimates of «i n f and M acc 
are ~ 2 km s" 1 and ~ 3 x 10~ 2 M Q yr" 1 for G31.41+0.31, and ~ 1.2 km s" 1 and ~ 9 x 10" 3 M Q yr" 1 for 
G24.78+0.08 Al. Such large accretion rates could weaken the effect of stellar winds and radiation pressure and 
allow further accretion on the star. The values of T rot and -/Vch 3 cn, derived by means of the RD method, for 
both G31.41+0.31 and the sum of cores Al and A2 (core A of Codella et al. \TW7l in G24.78+0.08 are in the 
range 132-164 K and 2-8 x 10 16 cm -2 . For G31.41+0.31, the most plausible explanation for the apparent toroidal 
morphology seen in the lower K transitions of CH3CN (12-11) is self-absorption, which is caused by the high 
optical depth and temperature gradient in the core. 
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1. Introduction 

The formation of massive stars represents a puzzle from a 
theoretical point of view. Unlike low-mass stars, they are 
believed to reach the zero-age main sequence still deeply 
embedded in their parental cores: in particular, Palla & 
Stahlcr (1993) predict that this occurs for stellar masses 
in excess of 8 Af . Once the star has ignited hydrogen 
burning, further accretion should be inhibited by radiation 
pressure and powerful stellar winds, with the consequence 

Send offprint requests to: M. T. Beltran, e-mail: 
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* Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by INSU/CNRS 
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that stars more massive than 8 M Q should not exist. Two 
formation scenarios have been proposed to solve this para- 
dox: non-spherical accretion (Nakano et al.Q.995 Jijina & 
Adams 1996) and merging of lower mass stars (Bonncll et 
al. Q998). Discriminating between these two models rep- 
resents a challenging observational goal. 

In this context, the detection of disks would strongly 
favour the accretion scenario, as naively one would ex- 
pect that random encounters between merging stars do 
not lead to axially symmetric structures since the disks 
associated with the lower mass stars should be destroyed 
during the merging process. On the contrary, conservation 
of angular momentum is bound to cause flattening and 
rotation of the infalling material, thus producing disk-like 
bodies. Indeed, circumstellar disks have been detected in 
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low-mass stars and found to undergo Keplerian rotation 
(Simon et al. 2000). Similar evidence has been found in a 
few high-mass Young Stellar Objects (YSOs), but in most 
cases the angular resolution was insufficient to assess the 
presence of a disk unambiguously. In conclusion, only few 
bona fide examples are known (see Cesaroni 2002 ) and all 
of these are associated with moderately massive stars (Bl- 
B0). This is not sufficient to understand the role of disks 
in the formation of even more massive stars and establish 
the relevance of accretion to this process. 

With this in mind, we have decided to perform a search 
for disks in high-mass YSOs by means of high angular res- 
olution observations. For this purpose, we have selected 
luminous (Lboi > 10 4 L Q ) objects with typical signposts 
of massive star formation such as water masers and ultra- 
compact (UC) Hn regions. The first objects selected for 
this project are G31.41+0.31 and G24.78+0.08. The for- 
mer is a well studied hot core located at 7.9 kpc (Cesaroni 
et al. 119941 I1998J1 . where preliminary evidence of a ro- 
tating massive disk oriented perpendicularly to a bipolar 
outflow has been reported in Cesaroni et al. I|1994|) and 
Olmi et al. (1996a). The latter one is a cluster of massive 
(proto)stars with a distance of 7.7 kpc, where recently 
Furuya et al. I|2002[> have detected a pair of cores, each 
of these associated with a compact bipolar outflow. On 
the basis of previous experience with this type of objects 
(see e.g. Cesaroni et al. [1999), CH 3 CN has been used as 
disk tracer. This is a low-abundance molecule which is ex- 
cited in very dense regions. Therefore, searching for disks 
requires not only high angular resolution, but also great 
sensitivity given the faintness of the lines observed. The 
first result derived from the analysis of the high angular 
resolution data has been the detection of clear velocity 
gradients perpendicular to the axis of the molecular out- 
flows in the cores (see Fig. 2 in Beltran et al. 2 0041 here- 
after Paper I), which points out the presence of rotating 
structures in the hot cores: one in G31.41+0.31 and three 
in G24.78+0.08. In Paper I, we referred to these rotat- 
ing structures as "disks"; however, their sizes (~ 8000- 
16000 AU in diameter) and masses (~ 80-500 Mq) make 
them remarkably different from the geometrically thin cir- 
cumstellar disks seen in low-mass YSOs, where the mass 
of the disk is much less than that of the star. Thus, these 
structures are more toroidal, and it is plausible to inter- 
pret such toroids as "circumcluster" disks, in the sense 
that they may host not just a single star, but a whole clus- 
ter. The rotating toroidal structures detected in Paper I 
are unstable as suggested by the fact that the mass of the 
cores is much larger than the dynamical mass needed for 
equilibrium. The main finding of the preliminary results 
of such a study is that nonspherical accretion is a viable 
mechanism to form high-mass stars. 

In this article we present a full report on the results 
obtained in Paper I, and a complete analysis of the data. 
We thoroughly analyze all the molecular lines detected 
toward the cores, the gas and dust distribution, and the 
physical parameters. We study in detail the kinematics of 
the cores and model the gas emission by assuming that 



it arises from a disk-like structure with an internal radial 
velocity field, and by adopting power-law distributions for 
the velocity and temperature of the emitting gas. 

2. Observations 

2.1. BIMA observations 

Millimeter observations of G24.78+0.08 were carried out 
with the BIMA array 1 in the A, B, C, and D configu- 
rations between 2002 September and 2003 February (see 
Table 2]). The digital correlator was configured to observe 
simultaneously the continuum emission and some molecu- 
lar lines. Details of the observations are given in Tabled 
The phase center was located at a(J2000) = 18 h 36 m 
12*659, <5(J2000) = -07° 12' 1C/15. The bandpass cali- 
bration was carried out on 3C345. Amplitude and phase 
were calibrated by observations of 1743—048, whose flux 
density of 4.4 Jy at 3.20 mm (A and B configurations) and 
of 4.7 Jy at 3.16 mm (C and D configurations) was deter- 
mined relative to Uranus, or Mars for the B configuration 
observations. The data were calibrated and analyzed us- 
ing standard procedures in the MIRIAD software package. 
We subtracted the continuum from the line emission di- 
rectly in the (u, ijj-domain in order to reduce non linearity 
effects in the deconvolution, and thus any amplification of 
errors induced in this process. 

2.2. PdBI observations 

Interferometric observations of G31.41+0.31 and 
G24. 78+0.08 were carried out with the IRAM Plateau 
de Bure Interferometer (PdBI) in the most extended 
configuration on 2003 March 16. By using the dual 
frequency capabilities of the PdBI we observed the 
continuum and molecular lines simultaneously at 3.3 and 
1.4 mm. The receivers were tuned single-sideband at 
3.3 mm and double-sideband at 1.4 mm. The correlator 
was centered at 92.0475 GHz in the lower sideband at 

3.3 mm and at 220.5250 GHz in the upper sideband at 

1.4 mm. The frequency setup of the correlator and the 
list of the observed molecular lines are shown in Tabic |5J 
The eight units of the correlator were placed in such a 
way that a frequency range free of lines could be used to 
measure the continuum flux. The phase center was set to 
the position a(J2000) = 18 h 47 m 34?330, c5(J2000) = -01° 
12' 46'/50 for G31.41+0.31, and to a(J2000) = 18 h 36 m 
1X661, <5(J2000) = -07° 12' 10'.'15 for G24.78+0.08. The 
bandpass of the receivers was calibrated by observations 
of the quasar 3C273. Amplitude and phase calibrations 
were achieved by monitoring 1741—038, whose flux 
densities of 4.63 Jy and 2.99 Jy at 3.3 mm and 1.4 mm, 
respectively, were determined relative to MWC349. The 
data were calibrated and analyzed with the GILDAS 
software package developed at IRAM and Observatoire 

1 The BIMA array is operated by the Berkeley-Illinois- 
Maryland Association with support from the National Science 
Foundation. 
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Synthesized beam 

Observation Configuration Frequency HPBW P.A. Bandwidth Spectral resolution rms noise a 

(GHz) (arcsec) (deg) (MHz) (km s -1 ) (mjy beam -1 ) 

3.20 mm continuum A, B 93.55 1.28 x 0.94 1200 - I 

CH 3 CN (5-4) A, B 91.98 0.83 x 0.60 61 25 0.31 50 

3.16 mm continuum C, D 94.75 19.8 x 16.3 -2 800 - 3 

N 2 H+ (1-0) C, D 93.18 14.3 x 13.0 16 25 0.31 100 

C 34 S (2-1) C, D 96.41 12.6 x 13.5 14 25 0.31 100 

CH 3 OH (2fc-lfc) G_D 96.74 14.0 x 12.6 12 25 031 100 

(a) For the molecular line observations the 1-a noise is per channel. 



Table 2. Frequency setups used for the molecular lines observed with the IRAM PdBI 



Line 


Center frequency 


Bandwidth 


Spectral resolution 


rms noise a 




(MHz) 


(MHz) 


(MHz) 


(km s" 1 ) 


(mjy beam -1 ) 


CH 3 CN (5-4), CH 3 l;i CN (5-4) 


91987.094 


80 


0.156 


0.5 


12 


CH 3 13 CN (5-4) 


91913.398 


320 


2.5 


8.1 


3 


CH3CN (5-4) vs = 1 


92353.438 


160 


0.625 


2.0 


7 


CH 3 CN (5-4) v a = 1 


92175.484 


320 


2.5 


8.1 


3 


CH 3 CN (12-11) 


220747.266 


80 


0.156 


0.2 


50 


CH3CN (12-11), CH 3 13 CN (12-11) 


220679.297 


80 


0.156 


0.2 


50 


CH 3 CN (12-11), CH 3 13 CN (12-11) 


220594.422 


160 


1.25 


1.7 


20 


CH 3 CN (12-11), CH 3 13 CN (12-11), 13 CO (2-1) 


220475.812 


320 


2.5 


3.4 


16 



(a) For the molecular line observations the 1-a noise is per channel. 



de Grenoble. We subtracted the continuum from the line 
emission directly in the (u, vj-domain. The synthesized 
CLEANed beams for maps made using natural weighting 
were 2'/7 x l'/2 at P.A. = -174° (3.3 mm) and l'.'l x 0'.'5 at 
P.A. = -170° (1.4 mm) for G31.41+0.31, and 3.0'.'l x l'/2 
at P.A. = -177° (3.3 mm) and l'/2 x 0"5 at P.A. = -174° 
(1.4 mm) for G24.78+0.08. Note that the G24.78+0.04 
BIMA and PdBI continuum data at 3 mm were merged in 
the {u, ?;)-domain (see next section) and the synthesized 
CLEANed beam for the resulting map was l'/3 x 0'.'9 at 
P.A. = 0°. 

One of the lines observed and clearly detected toward 
both regions is 13 CO (2-1). However, in both cases the 
line profile exhibits a strong lack of emission at the center, 
probably due to a missing flux problem since the interfer- 
ometer might have filtered out part of the extended emis- 
sion as suggested by the single-dish spectra (G31.41+0.31: 
Cesaroni, private communication; G24.78+0.08: Cesaroni 
et al. 12003(1 . In addition, although the line exhibits promi- 
nent wings, it has been impossible to disentangle the core 
emission (at a systemic velocity) from the high-velocity 
emission to derive the parameters of the outflowing mate- 
rial. Therefore, the emission of 13 CO (2-1) has not been 
analyzed in the paper. 

2.3. Merging the BIMA and PdBI data for 
G24.78+0.08 

The BIMA (A and B configuration) and PdBI contin- 
uum data of G24.78+0.08 at 3 mm were merged in the 
(u, u)-domain. The projected baseline lengths range from 
~ 3.8 to 416 kA for the BIMA data, and from ~ 4.8 to 



127 kA for the PdBI data. The BIMA (u, w)-data were 
first converted to FITS format using a procedure avail- 
able in the MIRIAD software package, and then converted 
to GILDAS format using a procedure available in the 
GILDAS software package. All the BIMA data were then 
merged together to create a single (u, u)-table using the 
task uv-merge of GILDAS, and the relative weights of the 
data were divided by 1 x 10 6 so that the weight density of 
the BIMA data matches that of the PdBI data. Finally, 
the PdBI (u, ?;)-data were merged with the already created 
BIMA (u, i;)-table. 

3. G31.41+0.31 

3.1. Results 

3.1.1. Continuum emission 

In Fig. ^ we show the PdBI maps of the 3.3 and 1.4 mm 
continuum emission toward the hot core G31.41+0.31. 
The continuum dust emission is marginally resolved at 
both wavelengths and shows a compact source plus 
an extended envelope. With the BIMA interferometer, 
Wyrowski et al. (2001) have mapped a similar structure 
at 1.4 mm. The position and fluxes at 3.3 and 1.4 mm, 
and the deconvolved size of the source, measured as the 
geometric mean of the major and minor axes of the 50% 
of the peak contour at 1.4 mm, are given in Table [3] The 
flux densities of 0.67 Jy at 3.3 mm and 4.4 Jy at 1.4 mm 
are consistent with the values of 0.28 Jy at 3.4 mm and of 
4.0 Jy at 1.3 mm measured by Maxia et al. I|2001|l with the 
Owens Valley Radio Observatory (OVRO) interferometer. 
Hatchell et al. ( 2000| have measured a flux of 4.9 Jy at 
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Table 3. Positions, flux densities and sizes at millimeter wavelengths for the cores in G31.41+0.31 and G24.78+0.08 





Position a 


/r ak (3.3mm) 
(Jy beam ) 




JP cak ( 1.4mm) 
(Jy beam ) 




Source Diameter b 


Core 


a(J2000) 
h m s 


5(J2000) 

o / // 


S„(3.3mm) 

(Jy) 


S„ (1.4mm) 

(Jy) 


e 3 

(arcsec) 


0s 

(AU) 


G31 


18 47 34.32 


-01 12 45.9 


0.252 


0.674 


1.3 


4.4 


1.1 


8700 


G24 Al 


18 36 12.57 


-07 12 10.9 


0.086 c 


0.211 d 


0.191 


0.882 


1.2 


9200 


G24 A2 


18 36 12.49 


-07 12 10.0 


0.086 c 


0.211 d 


0.125 


0.574 


0.9 


7000 


G24 B 


18 36 12.66 6 


-07 12 15. 2 C 


0.016 


0.040 


< 0.020 f 


< 0.020 f 


1.5 s 


11500 E 


G24 C 


18 36 13.12 


-07 12 07.4 


0.005 


0.005 


0.056 


0.100 


0.9 


7000 



(a) Position of the 1.4 mm emission peak. 

(b) Deconvolved geometric mean of the major and minor axes of the 50% of the peak contour at 1.4 mm. 

(c) Peak intensity of core A (cores Al and A2 together; see Sect. 14. Ol , 

(d) Flux density of core A (cores Al and A2 together; see Sect. 14. Oil . 

(e) Source not detected at 1.4 mm. Position corresponds to the 3.3 mm emission peak. 

(f) Upper limit for nondetected source is 3er. 

(g) Deconvolved geometric mean of the major and minor axes of the 50% of the peak contour at 3.3 mm. 

G31.41+0.31 continuum 



-1°12'42" 
-1°12'44" 
g" -1°12'46" 

o 

C\l 

<o -1°12'48" 



-1°12'50" 
-1°12'52" 

18 h 47 m 34 s 6 34 s 4 34 s 2 34 s 18 h 47 m 34 s 6 34 5 4 34 s 2 34 s 

a(J2000) a(J2000) 

Fig. 1. PdBI map of the 3.3 mm (left panel) and the 1.4 mm (right panel) continuum emission toward the core of 
G31.41+0.31. The contour levels range from 20 to 240 mJy beam -1 in steps of 20 mJy beam -1 for the 3.3 mm map, 
and from 100 to 1300 mJy beam -1 in steps of 200 mJy beam -1 for the 1.4 mm map. The synthesized beam is shown 
in the lower left-hand corner. The white dots mark the position of the 7 mm continuum emission peaks detected by 
Homer (private communication). 




1.35 mm inside the 22" James Clerk Maxwell Telescope 
(JCMT) beam. The fact that the continuum flux seen with 
a single dish is not much higher than that measured with 
the interferometers seems to indicate that the source is 
very centrally peaked, and that there is no significant frac- 
tion of the total flux filtered out by the interferometers. 
In Paper I we derived a core mass of 490 Af© from the 
millimeter continuum emission assuming a mass opacity 
of ~ 0.02 cm -2 g -1 at 1.4 mm for a gas-to-dust ratio of 
100 and a temperature of 230 K (Olmi ct al. I199fiajl . 

3.1.2. CH 3 CN and CH 3 13 CN 

Figure HI shows the spectra of the CH 3 CN and CH 3 13 CN 
(12-11) and (5-4) lines toward G31.41+0.31. The first vi- 
brational state above the ground of the CH 3 CN (5-4) line, 
which is denoted Vg = 1, is also shown in this figure. The 
spectra have been taken integrating the emission over the 
3cr contour level area. Figure shows the maps of the 



CH3CN (12-11) emission averaged under the K = 1 and 
2 components and under the K = 8 component, and the 
CH3 13 CN (12-11) emission averaged under the K = 1 and 
2 components. Figure 0] shows the maps of the CH3CN 
(5-4) and the CH3 13 CN (5-4) emission averaged under 
the K — to 4 components, and the CH 3 CN (5-4) v s = 1 
emission averaged under all observed transitions. 

Several X-components of the different rotational tran- 
sitions of methyl cyanide are clearly detected toward 
G31.41+0.31 as can be seen in Fig. □ We fitted multi- 
ple i^-components of a given rotational transition by fix- 
ing their separation in frequency to the laboratory values 
(see Lovas l*H)92l and references therein) and forcing their 
line width to be identical. As different if-components of 
a same rotational transition fall in different PdBI correla- 
tor units with different spectral resolutions, sometimes we 
did not fit simultaneously all of the lines for a given ro- 
tational transition, but groups of them. The groups were 
of 5 lines at maximum because CLASS, which is the soft- 
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ware package used to analyze the lines, only allows to fit 
5 lines simultaneously. In addition, in a few cases, the 
line width in the fit was fixed to the value obtained for 
other ^-components of the same transition; otherwise, 
due to the low spectral resolution with which those tran- 
sitions have been observed, the line widths obtained when 
fitting the profiles were unrealistically wide. The fact of 
fixing the line width for groups of lines with different 
if-components fitted simultaneously does not affect sig- 
nificantly the results and conclusions derived from the 
analysis of the lines. The profiles show no evidence of 
self-absorption, are quite symmetric and can be reason- 
ably fitted by a Gaussian profile. The line parameters of 
the ground state transitions of CH 3 CN and CH 3 13 CN are 
given in Table 01 while those of the v& = 1 are given in 
Table Note that in some cases the lines are so blended 
with other methyl cyanide components (e.g. the CH3 13 CN 
K = with CH 3 CN K = 5 for the (12-11) transition), 
or with other molecular lines such as 13 CO (2-1), that it 
has been impossible to derive their parameters. As can 
be seen in Tables 0] and [5] there is a good agreement be- 
tween the systemic velocities derived from all transitions 
(Vlsr — 97 km s _1 ). Regarding the line widths, we found 
a similar value, FWHM ~ 8 km s _ \ for CH 3 CN (12-11), 
CH 3 13 CN (12-11), and CH 3 CN (5-4) v s = 1. On the other 
hand, for CH 3 CN (5-4) and CH 3 13 CN (5-4) we found val- 
ues of ~ 10 and ~ 6 km s , respectively. Assuming that 
the line profiles are Gaussian, the ratio of the observed 
line width (broadened by the line optical depth) and the 
intrinsic one, can be written as 



A^obs 

AV5„t 



1 



-In 



In 



1 + c-- T ° 



(1) 



where tq is the optical depth at the center of the line. In 
practice, Eq. {JJ gives the line width ratio between an op- 
tically thick and an optically thin transition of the same 
species, where to corresponds to the optical depth of the 
thick line. One may hence apply Eq. (JIJ to the K lines of 
CH 3 CN (5-4) and CH 3 13 CN (5-4). The optical depth r 
can be derived from the brightness temperature ratio be- 
tween lines of the two species with the same K (assuming a 
relative abund ance [CH 3 CN]/[CH 3 13 CN]=50; see Wilson 
& Rood 199|J: using this in Eq. one obtains the line 
width ratio between CH 3 CN and CH 3 13 CN expected if 
line broadening is due only to optical depth effects. Such 
a ratio turns out to range between 1.8 and 2.1, roughly 
consistent with that measured directly from the line pro- 
files: 1.6 ± 0.02. We conclude that the CH 3 CN lines are 
wider than the CH 3 13 CN ones as a consequence of the 
larger optical depth, although one cannot exclude that 
different methyl cyanide transitions are tracing gas with 
slightly different physical conditions as well. 

As can be seen in Fig. 1 of Paper I, and in Fig. the 
CH 3 CN K = 1 and 2 components seem to trace a toroidal 
structure with the millimeter continuum emission peak- 
ing at the central dip. To explain this morphology two 
possible scenarios are suggested: (i) a dramatic CH 3 CN 




220.3 220.4 220.5 220.6 220.7 
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00 1 




I 1 


[ 1 1 
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1 , 4 


i , i \i 
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91.9 91.92 91.94 91.96 91.98 



0.5 




92.15 



92.2 



92.25 

v (GHz) 



92.3 



92.35 



Fig. 2. Methyl cyanide spectra integrated over the 3cr con- 
tour level area toward G31.41+0.31 as seen with the PdBI. 
We show in the top CH 3 CN (12-11), in the middle CH 3 CN 
(5-4), and in the bottom CH 3 CN (5-4) vibrationally ex- 
cited (i>8 = 1). Different K numbers (top and middle pan- 
els) are marked with dashed lines in the upper (lower) part 
of each spectra in the case of CH 3 CN (CH 3 13 CN). Note 
that different if-components of a same transition have dif- 
ferent spectral resolution because of the different spectral 
resolution of the correlator units (see Table 0). 



abundance drop in the central region of the core; or (ii) 
self-absorption, due to the high optical depth at the center 
of the core and the existence of a temperature gradient. 
Figures and [5] show that the different if-components 
do not trace the same material, and that for large val- 
ues of K the CH 3 CN emission peaks closer to the center. 
In fact, as shown in Fig. [5J the K = 8 emission fills the 
central dip visible in the K = 0, 1, and 2 emission. The 
CH 3 13 CN K —1 and 2 emission is also found closer to 
the center than the CH 3 CN K =1 and 2 emission. Thus, 
the lack of emission at the center of the core occurs only 
for the lower if -components of CH 3 CN, unlike the higher 
if-components, CH 3 13 CN, and continuum emission. This 
indicates that the CH 3 CN molecules are not destroyed at 
the center of the core. As can be seen in Fig-El the CH 3 CN 
(12-11) K = to 4 lines show similar intensities, sugges- 
tive of high optical depth. Note that \aw-K transitions 
of CH 3 13 CN, which have excitation energies very similar 
to those of the corresponding CH 3 CN K lines, trace a 
slightly inner region: this is an opacity effect as the less 
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Table 4. CH 3 CN and CH 3 13 CN line parameters for G31.41+0.31 



Line 


Vlsr 


FWHM 


n 


j T B dV 


rms 


K 


(km s" 1 ) 


(km s" 1 ) 


(K) 


(K km s- 1 ) 


(K km s" 1 ) 


CH 3 CN (5-4) 





96.29 ±0.03 


9.77 ±0.04 


19.0 ±1.3 


198 ±3 


0.8 


1 


33 


3? 


15.9 ±1.3 


165 ±3 




2 


33 


33 


21.6 ± 1.3 


225 ±2 




3 


33 


33 


22.2 ± 1.3 


231 ± 2 


33 


4 


33 


33 


14.2 ± 1.3 


148 ±2 




CH 3 ia CN (5-4) 





97.41 ± 0.03 


6.11 ±0.05 


5.83 ±0.46 


37.9 ± 0.6 


0.24 


1 


33 


33 


6.45 ± 0.46 


41.9 ± 0.6 


33 


2 


n 


33 


4.76 ± 0.46 


31.0 ±0.6 


33 


3 


35 


33 


5.41 ±0.46 


35.2 ±0.6 


33 


4 


97.41 b 


6.11 b 


2.17±0.18 c 


14.1 ± 1.2 d 


0.08 


CH3CN (12-11) 





96.24 ±0.06 


8.07 ±0.08 


36.7 ± 1.4 


315 ±8 


2.9 


1 


33 


33 


31.3 ± 1.4 


269 ±8 


33 


2 


33 


33 


39.0 ± 1.4 


335 ±8 


33 


3 


33 


33 


40.4 ± 1.4 


347 ±8 


33 


4 


96.84 ± 0.09 


8.10 ±0.12 


34.9 ± 1.2 


301 ±7 


2.5 


5 


33 


33 


26.9 ± 1.2 


232 ± 11 


33 


6 


96.78 ±0.16 


8.00 b 


24.3 ±3.0 


207 ±7 


1.5 


7 


97.15 ±0.35 


8.00 b 


16.1 ±3.0 


137 ± 16 


3.6 


8 


33 


33 


10.9 ±3.0 


93 ± 18 


33 


9 


97.85 ±1.28 


8.00 b 


(e) 


(e) 


3.4 


10 


33 


33 


3.4 ±0.7 


29 ± 22 


33 


11 


97.12 ±0.65 


8.00 b 


1.0 ± 1.6 


8 ±17 


2.6 


CH 3 l;i CN (12-11) 





96.84 ±0.09 


8.10 ±0.12 


(f) 


(f) 


2.5 


1 


33 


33 


14.2 ± 2.5 


122 ±8 


33 


2 


33 


33 


12.9 ±2.5 


111 ± 6 


33 


3 


96.78 ±0.16 


8.00 b 


19.8 ±3.0 


169 ±7 


1.5 


4 


33 


33 


9.0 ±3.0 


76.3 ±6.6 


33 


5 


33 


33 


5.0 ± 3.0 g 


42.2 ± 7.2 g 


33 


6 


33 


33 


3.8 ±3.0 


32.1 ± 6.6 


33 


7 


97.85 ±1.28 


8.00 b 


2.6 ±0.7 


22 ±21 


3.4 


8 


97.12 ±0.65 


8.00 b 


15.0 ± 1.6 h 


127 ± 17 h 


2.6 


9 


33 


33 


6.4 ± 1.6 


55 ± 17 


33 



(a) Brightness temperature integrated over a 3a emission level area. 

(b) This parameter was held fixed in the fit. 

(c) Computed from the integrated line intensity. 

(d) Measured with the command print area of CLASS. 

(e) Blended with 13 CO (2-1). 

(f) Blended with CH3CN K = 5. 

(g) Blended with CH 3 CN K = 7. 

(h) Dominated by wing emission of 13 CO (2-1). 



abundant isotopomer traces higher column densities. In 
order to confirm a possible opacity effect, we computed 
the optical depth at each point using the K — 2 compo- 
nents of CH 3 CN and CH 3 13 CN, assuming that CH 3 CN is 
optically thick and that the excitation temperature is the 
same for both transitions. A map of the optical depth of 
the CH3CN main species, T12, shows that the opacity in- 
creases toward the center of the core: it ranges from 10 to 
50 for K ~ 2, and from 5 to 20 for K = 6. In conclusion, 
the most plausible explanation for the apparent toroidal 
geometry of the CH3CN map in Fig. is self-absorption, 



which is caused by the high optical depth and the tem- 
perature gradient in the core. In Sect. I3~2~T1 we discuss in 
more detail this temperature gradient. 

3.2. Discussion 

3.2.1. Temperature gradient 

In Sect. I3.L21 we have illustrated the existence of a tem- 
perature gradient in the G31.41±0.31 core. Here we want 
to estimate the temperature and its variation across the 
core. Following Olmi et al. Q1993P we derived estimates 
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Table 5. CH 3 CN (5-4) v 8 = 1 line parameters for G31.41+0.31 



Line 


Vlsr 
(km s 1 ) 


FWHM 
(km s 


n 

(K) 


J'T B dV 
(K km s" 1 ) 


rms 
(K km s _1 ) 


(1,1) 


97.45 ±0.41 


8.02 ± 0.44 


5.15 ±0.24 


44 ±5 


1.0 


(2,1) 


)} 


11 


4.69 ±0.24 


40 ±5 


11 


(0,1) 


11 


11 


4.62 ±0.24 


39 ±5 


11 


(3,1) 




11 


3.65 ±0.24 


31 ±5 


11 


(1,-1) 




11 


4.12 ±0.24 


35 ±5 


11 


(4,1) 


97.45 ± 0.50 


8.71 ± 0.69 


3.17 ±0.31 


29 ±7 


1.2 


(2,-1) 


)} 


11 


4.95 ±0.31 


46 ±7 


11 


(3,-1) 


11 


11 


3.27 ±0.31 


30 ±7 


11 


(4,-1) 


11 


11 


1.63 ±0.31 


15 ±6 


11 


(1,1) 


97.45 b 


8.71 b 


4.72±0.12 c 


39 ± l d 


0.09 



(a) Brightness temperature integrated over a 3a emission level area. 

(b) This parameter was held fixed in the fit. 

(c) Computed from the integrated line intensity. 

(d) Measured with the command print area of CLASS. 



for the rotational temperature, T rot , and the total methyl 
cyanide column density -/Vch 3 cn, by means of the ro- 
tation diagram (RD) method, which assumes that the 
molecular levels are populated according to LTE condi- 
tions at a single temperature T rot . In the high density limit 
where level populations are thermalised, one expects that 
T m t — Ik) the kinetic temperature. The CH3CN ground 
level transitions appear to be optically thick, as suggested 
by the ratio between main species and isotopomer. As 
already mentioned in the previous section, T12, which is 
position depending and has been evaluated at the sys- 
temic velocity, ranges from 10 to 50 for K = 2, and from 
5 to 20 for K — 6 in the core. Thus, in the fit per- 
formed to the Boltzmann plot we used only the transi- 
tions of CH 3 13 CN, and those from the CH 3 CN v s = 1 
state (see Fig. EJ). All the spectra were integrated over a 
region of ~ 4" 7 in diameter (~ 17 arcsec 2 ), in order to 
include all the emission inside the 3er contour level for 
the (5-4) transition. Assuming that also the emission due 
to the vibrational state is optically thin and a relative 
abun dance [CH 3 CN]/[CH 3 13 CN]=50 (see Wilson & Rood 
1994), the T rot and source averaged -/Vch 3 cn derived by 
taking into account all transitions mentioned above are 
164 ± 5 K and 8 x 10 16 cm -2 , respectively. We also es- 
timated T ro t and Ncn a cN from each set of J ± 1 — > J 
lines separately. The rotation temperatures estimates are 
135 ± 2 and 220 ± 16 K, and the column densities 9 x 10 16 
and 7 x 10 16 cm~ 2 for the (5-4) and (12-11) transitions, 
respectively. These values of T rot are consistent with those 
derived from CH 3 CN and CH 3 13 CN (6-5) by Olmi et 
al. Ijl996a|l for the compact component of the molecular 
clump (~ 1-2"), called "core", and the most extended one 
(~ 10-20"), called "halo", which are 230 and 130 K, re- 
spectively. The fact that the T rot that we derived by using 
only the (12-11) transitions is close to the value derived by 
Olmi et al. (19963) f° r the "core", and that derived using 
only the (5-4) lines is closer to the value for the "halo" , 
could be indicating that the highest excitation transitions 
are tracing material closer to the embedded source, while 



the lowest excitation transitions are tracing material in 
the outer shell of the core. From this analysis, it is clear 
that the temperatures in G31.41±0.31 are very high, con- 
firming that the methyl cyanide emission comes from hot 
molecular gas. Regarding the column densities, the esti- 
mates of the order of ~ 10 17 cm -2 are consistent with 
those obtained by Olmi et al. l|1996al[T996bll . 

In order to establish the existence of a temperature 
gradient in the core, as suggested by Figs. and \5\ (see 
Sect. 13.1.21) . we computed a map of T rot by fitting the 
Boltzman plots at each position (pixel) in the core. The 
fits were performed with the RD method to the CH 3 13 CN 
K=l, 2, 3, 4, 5, and 6 components (K=0 was not used be- 
cause too blended with CH 3 CN K=5). Figured shows the 
resulting map of T rot , as well as the total column density 
map, -/Vch 3 cn, derived with the same method; for the sake 
of comparison the CH 3 13 CN (12-11) emission averaged 
under the K = 1 and 2 components is also shown. This fig- 
ure clearly confirms that the temperatures reached toward 
the center of the core, where the millimeter continuum 
source is embedded and the lower if-components methyl 
cyanide emission shows a dip, are very high (~ 350 K), 
while the temperatures toward the peak of CH 3 13 CN are 
lower (~ 200-250 K). The temperature at the position 
of the peaks of the torus seen in the CH 3 CN emission 
(see Fig. 01 is even lower, T rot ~ 150 K. For this source 
Cesaroni et al. (|1998|) have derived a temperature gradient 
of the type T oc i? -3 / 4 using NH 3 (4,4): this is consistent 
with our findings 350 ^ 150 x 3 3 / 4 , where 3 is roughly 
the value of the distance where the temperature is 150 K 
divided by the distance where the temperature is 350 K. 
The ATch 3 cn map also shows that the gas column density 
increases toward the center, with the peak at the position 
of the millimeter continuum source, confirming that the 
CH 3 CN gas is present also in the central region of the 
core. 

The geometry of the G31.41±0.31 core seems to be 
the result of the combination of two effects: the exis- 
tence of a temperature gradient in the core plus opac- 
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Fig. 3. PdBI intensity maps of the CH 3 CN (12-11) emis- 
sion averaged under the K = 1 and 2 components (top 
panel), under the K — 8 component (middle panel), 
and of the CH3 13 CN (12-11) emission averaged un- 
der the K = 1 and 2 components (bottom panel) to- 
ward G31.41+0.31. The contour levels range from 0.1 to 
0.94 Jy beam -1 (top panel), to 0.70 Jy beam -1 (middle 
panel), and to 0.82 Jy beam -1 (bottom panel), in steps 
of 0.12 Jy beam -1 . The synthesized beam is shown in 
the lower left-hand corner. The white cross marks the po- 
sition of the 1.4 mm continuum emission peak, and the 
white dots the 7 mm continuum emission peaks detected 
by Hofner (private communication). 
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Fig. 4. PdBI intensity maps of the CH 3 CN (5-4) (top 
panel) and the CH3 13 CN (middle panel) emission aver- 
aged under the K = to 4 components, and of the vg = 1 
emission (bottom panel) averaged under the (1,1), (2,1), 
(0,1), (3,1), (1,-1), (4,1), (2,-1), (3,-1), (4,-1), and (1,1) 
transitions toward G31.41+0.31. The contour levels range 
from 0.1 to 0.70 Jy beam -1 in steps of 0.12 Jy beam -1 (top 
panel), from 0.02 to 0.22 Jy beam -1 (middle panel), and to 
0.18 Jy beam -1 (bottom panel) in steps of 0.04 Jybeam -1 . 
The synthesized beam is shown in the lower left-hand cor- 
ner. The white cross marks the position of the 1.4 mm con- 
tinuum emission peak, and the white dots the 7 mm con- 
tinuum emission peaks detected by Hofner (private com- 
munication) . 



ity. In Sect. 13. Ol we checked for a possible opacity effect 
by computing the optical depth at each point assuming, 
as a first approximation, a single excitation temperature. 
However, because of the existence of a temperature gradi- 
ent in the core, the optical depth can be better measured 
taking it into account. In order to better fit the data and 
understand the geometry of the core we computed models 



which solve the radiative transfer equation. These models 
are discussed in the next section. 

3.2.2. Velocity field: the rotating toroid 

In Paper I we detected a clear velocity gradient in 
the G31.41+0.31 core by simultaneously fitting multiple 
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G31 .41 +0.31 




a (J2000) 

Fig. 5. Overlay of the PdBI CH 3 CN (12-11) emission 
averaged under the K = 0, 1, and 2 components (con- 
tours) on the emission under K — 8 (greyscale) to- 
ward G31.41+0.31. Contour levels range from 0.1 to 
0.94 Jybeam -1 in steps of 0.12 Jybeam -1 . Greyscale 
levels range from 0.1 to 0.70 Jybeam -1 in steps of 
0.12 Jybeam -1 km s . The synthesized beam is shown 
in the lower left-hand corner. The white cross marks the 
position of the 1.4 mm continuum emission peak, and the 
white dots the 7 mm continuum emission peaks detected 
by Hofner (private communication). 
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Fig. 6. Rotation diagram for G31.41+0.31 with superim- 
posed fit. Filled circles, open circles, and filled triangles 
correspond to the CH 3 13 CN (5-4), the CH 3 13 CN (12-11), 
and the CH3CN (5-4) Ug = 1 transition, respectively. 

CH3CN (12-11) i^-components at each position where 
CH3CN is detected. This velocity gradient is perpendic- 
ular to the molecular outflow detected in the region by 
Olmi et al. (1996a), and as already discussed in Paper I, 
the most plausible explanation for such a velocity gradient 
is that the core is undergoing rotation about the outflow 
axis. Also Gibb et al. I|2004fl have studied the velocity 
field in G31.41+0.31 by analyzing the C 18 and H 2 S line 
emission. By mapping the H 2 S line wings, these authors 
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G31.41 - CH 3 13 CN (12-11) 




a(J2000) 

Fig. 7. Overlay of the PdBI CH 3 13 CN (12-11) emission 
averaged under the K = 1 and 2 components (contours) 
on the T rot map (top panel) and the iVcHsCN map (bot- 
tom panel) in greyscale, derived by fitting the Boltzmann 
plots at each position, toward G31.41+0.31. The con- 
tour levels range from 0.1 to 0.82 Jy beam -1 in steps of 
0.12 Jy beam -1 . Greyscale levels range from 100 to 350 K 
by 50 K (top panel), and from 2 to 7 x 10 17 cm -2 by 
1 x 10 17 cm -2 (bottom panel). The white cross marks the 
position of the 1.4 mm continuum emission peak, and the 
white dots the 7 mm continuum emission peaks detected 
by Hofner (private communication). 

have detected a velocity gradient in the same direction 
as our CH3CN gradient, but have interpreted this as a 
bipolar outflow, with the blueshifted emission toward the 
west and the redshifted toward the east. To our knowl- 
edge a high-density tracer such as CH3CN has never been 
found to trace outflows: instead there are clear examples of 
CH 3 CN emission tracing rotating disks, for example IRAS 
20126+4104 (Cesaroni et al . 113331 or G24.78+0.08 (Paper 
I); in the latter the velocity gradient seen in CH3CN in 
core Al is clearly perpendicular to the bipolar outflow 
mapped in 12 CO (Furuya et al l2002[) . We hence believe 
that the interpretation given in Paper I for the CH3CN 
emission in G31.41+0.31 is to be preferred to that of Gibb 
et al. (|2004ll . and consider the velocity gradient detected 
a signature of rotation. Only high-angular resolution ob- 
servations in unambiguous outflow tracers such as 12 CO 
may definitely prove which interpretation is correct. 

The velocity shift measured over an extent of ~ 
17000 AU is 4.2 km s -1 (Paper I), which corresponds to a 
velocity gradient of about 50 km s -1 pc -1 . The size and 
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the mass of ~ 500 M Q suggest that this toroid could be 
hosting not a single massive YSO but a cluster of YSOs. 
This is also suggested by the continuum emission at 7 mm 
toward the center, which reveals two point-like sources 
oriented in the NE-SW direction (Hofner, private com- 
munication), i.e. just along the toroid plane. It is hence 
tempting to speculate that collapse in the toroid has led 
to their formation. 

The kinematics of the gas toward the core of 
G31.41+0.31 can be seen in the position-velocity (PV) 
cut of the CH3CN (12-11) K = 3 emission along a di- 
rection (P. A. ~ 70°) where the velocity gradient is max- 
imum (Fig. IBJ. Clearly, Keplerian rotation, which has 
been detected toward lower luminous objects such as 
IRAS 20126+4104 (Cesaroni et al. 1997, 1999 Zhang et 
al. llTO), N GC 7538S (Sandell et aL lOTty . or M17 (Chini 
et al. 120040 , is not possible as the mass of the toroid is 
much larger than any reasonable stellar mass. In fact, the 
PV cut shows no hint of Keplerian rotation. However, it 
is interesting to see that according to theoretical models 
(Galli, private communication), discriminating between a 
velocity field with constant rotation velocity or with con- 
stant angular velocity should help to assess whether mag- 
netic braking, and thus magnetic field, plays a crucial role 
during the process of disk, and hence star, formation. In 
the absence of magnetic field, a self-gravitating isothermal 
quasistatic equilibrium situation is reached when the ro- 
tation velocity is constant, v — constant (7 = 0). In case 
that magnetic field is present, it couples the core with 
the surrounding medium by means of magnetic braking. 
This coupling yields a constant angular velocity, v oc R 
(7 = 1), in the whole core, which rotates as a rigid body 
and leads by ambipolar diffusion to an equilibrium situa- 
tion. Thus, in order to discriminate between the two possi- 
ble scenarios, which could give an indication of the stabil- 
ity of the cores and a hint on their evolutionary stage, we 
modeled the CH3CN emission by assuming that the emis- 
sion arises from a disk-like structure seen edge-on by the 
observer with an internal radial velocity field (v oc -R 7 ); 
a constant rotation velocity field and a constant angular 
velocity field. 

The models were computed by adopting a power law 
dependence on the distance from the core center also for 
the temperature T oc R~ q of the emitting gas. According 
to accretion disk theory for geometrically thin disks the 
dependence of the temperature on the radius is T oc i? -3 / 4 
(see Natta 120001 and references therein). Although geo- 
metrically thick, these toroids resemble more 2-D struc- 
tures like circumstellar disks than 3-D spherically sym- 
metric cores. Therefore, we have adopted a temperature 
profile T oc i?~ 3 / 4 in our models. This is also indicated 
by the findings of Cesaroni et al. (1998), who have prop- 
erly modeled the NH 3 (4,4) emission toward G31.41+0.31 
with T oc i? -3 / 4 . Regarding the density distribution, for 
simplicity we adopted constant density in the core. The 
synthetic PV diagrams were computed along the projected 
major axis of the toroid, and convolved with a Gaussian in 
the PV plane to obtain the same angular and spectral reso- 



Fig.8. Overlay of the PV plot of the CH3CN (12-11) 
K = 3 emission (grayscale) along the major axis of the 
toroid (P.A. ~ 70°) in G31.41+0.31, and the synthetic 
emission PV plot of a model (contours) with constant an- 
gular velocity (top panel), and with constant rotation ve- 
locity (bottom panel). The 2/-axis plots the difference be- 
tween the measured velocity and the systemic velocity, 
which is 97 km s -1 for G31.41+0.31. The contour levels 
range from 6 to 54 Jy beam -1 in steps of 8 Jy beam -1 . 
The error bars in the lower left-hand corner indicate the 
angular and spectral resolution of the data and the mod- 
els. 

lutions as those of the data. The parameters of the models 
are the inner and outer radius of the toroid, Ri nn and i? ou t ? 
respectively, a parameter related to the column density at 
the peak of the emission, the observed line width, AV, the 
power-law index p of the density, the temperature at the 
outer radius, T ou t, the power-law index q of the tempera- 
ture, and the rotation velocity measured at the edge of the 
toroid, v I0 t- The initial guesses for the input parameters 
such as the size of the toroid or the rotation velocity were 
derived from the observations (see Paper I) . 

The observed and synthetic PV diagrams for both sce- 
narios, constant angular velocity (angular) and constant 
rotation velocity (linear), can be seen in Fig. [S] From the 
modeling itself it is not possible to distinguish between 
constant angular (v oc R; top panel) or constant rota- 
tion velocity (v = constant; bottom panel), because both 
velocity fields suitably fit the data. The models were com- 
puted with and without temperature gradients. However, 
only those models with T oc i? -3 / 4 can fit properly the 
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Table 6. Best fit model parameters for G31.41+0.31 and core Al in G24.78+0.08 
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(a) Column density at the peak of the CH3CN emission. 
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Fig. 9. Line width plot of CH 3 13 CN (12-11) toward 
G31.41+0.31. The contour levels range from 3 to 7 km s _1 
in steps of 1 km s _1 . The white cross marks the position of 
the 1.4 mm continuum emission peak, and the white dots 
the 7 mm continuum emission peaks detected by Hofncr 
(private communication) . 



observed PV plots and reproduce the "inner hole" seen in 
the data. A similar result could be also obtained with a 
bigger Ri nn . However, in this case the hole at the center 
of the toroid would be also visible in the millimeter con- 
tinuum emission maps. The best fit model parameters are 
given in Table El The value of i? inn is 0/.T7 (1340 AU), 
which is roughly the separation between the two sources 
detected at 7 mm (Hofner, private communication). The 
values of v ro t and R ou t are consistent with those derived 
directly from the velocity gradient and the 1.4 mm contin- 
uum emission at 50% of the peak, respectively (see Table 1 
in Paper I). The value of the column density at the peak 
of the CH 3 CN emission, Aq^^j^, is ~ 20 times lower than 
the column density derived by means of the RD method 
(see previous section). This is due to the fact that the mod- 
els do not take into account the dumpiness of the region 
(i.e. a beam filling factor effect). The models also support 
the scenario according to which the lack of emission at 
the center of the toroid is not a geometric effect but is 
due to self-absorption, which is caused by the high tem- 
peratures and optical depths reached toward the center of 
the core, as suggested by the CH 3 CN and CH 3 13 CN inte- 
grated emission plots and by the column density distribu- 
tion obtained from the Boltzman plots (see Sect. l3.l!2l and 



13.2. 1|) . In fact, optically thin models with a temperature 
gradient cannot fit the data because in such a case, the 
external layers of the toroid become too optically thin to 
absorb enough photons from the center to reproduce the 
lack of emission or self-absorption visible in CH3CN. Thus, 
the models indicate that the emission in G31.41+0.31 has 
to be optically thick and that a temperature gradient is 
present in the core. 

The dynamical mass needed for equilibrium, Md yn , 
has been derived from the expression Afdyn = 
Vf 0t R ou t/G sin 2 i, where i is the inclination angle assumed 
to be 90° for an edge-on toroid. Note that this expression 
in Paper I is wrong: sin 2 i should be dividing, but the val- 
ues of Mdyn given in Table 1 were correctly computed. 
The value obtained is Md yn ~ 44 M Q . The fact that Md yn 
is much smaller than the mass of the core suggests, as al- 
ready pointed out in Paper I, that such structure is unsta- 
ble against gravitational collapse. Note that the magnetic 
field required to stabilize it would be of a few 20-40 mG, 
values too large to be plausible even in regions as dense 
as 10 8 cm" 3 (see Fig. 1 of Crutcher ITMIJl . 

In case that the toroid is undergoing collapse the gas 
should be more perturbed toward the center, where the 
1.4 mm continuum source is embedded, and the veloc- 
ity dispersion should increase inward. Thus, one would 
expect the line width, AV, to increase toward the cen- 
ter of the toroid. Indeed Fig. |5] shows that AU, which 
has been obtained by fitting the optically thin CH3 13 CN 
(12-11) line profiles at each position with a Gaussian, in- 
creases toward the center of the core, where the YSO(s) 
is embedded. One can obtain an estimate of the infall ve- 
locity, w; n f, assuming that the only contribution to the 
variation of AU is infall and that the contribution of tur- 
bulence is negligible. In such a case, i>; n f is given by the 
difference between AU at the position of the 1.4 mm 
continuum source and that at the edge of the toroid, 
V[ n f ~ (AUnn — AU out )/2. Interestingly, the value ob- 
tained is Ui n f ~ 2 km s , consistent with the i> ro t ob- 
tained from the models. Such agreement between u; n f and 
v ro t has been predicted by some recent theoretical stud- 
ies (see e.g. Allen et al. 2003:). Following the expression 
given in Paper I we estimated the accretion rate and ob- 
tained Af acc ~ 3 x 10~ 2 Mq yr _1 . Accretion rates that 
large have been estimated for the massive Class ob- 
ject IRAS 23385+6053 (Molinari et al. lTMg|) . and for the 
parsec-scale clumps in which high-mass star formation is 
observed (Fontani et al. l20tT2~fl , and could support the the- 
ories that predict that high-mass stars might form through 
non-spherical accretion with large accretion rates. 
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Fig. 10. Overlay of the BIMA integrated intensity maps of the N 2 H+ (1-0) emission (contours; top panel), the 
C 34 S (2-1) emission (contours; middle panel), and the CH3OH (2fc— lfc) emission (contours; bottom panel) on the 
BIMA 3.16 mm continuum emission (greyscale) toward G24.78+0.08. The contours are 3, 9, 15, 25, 45, 75, 120, and 
150 Jy beam- 1 km s" 1 for N 2 H+, 1.2, 3.6, 6, 10, 18, 30, and 48 Jy beam" 1 km s" 1 for C 34 S, and 3, 6, 12, 24, 48, 96, 
and 144 Jy beam" 1 km s" 1 for CH 3 OH (2 fe -l fe ), and 3, 9, 15, 25, 45, 60, and 120 times 2.5 mJy beam" 1 for the 
3.16 mm continuum emission. The dashed circumference in the top panel represents the BIMA primary beam (50% 
attenuation level). The synthesized beams are shown in the lower right-hand corner. The crosses pinpoint the positions 
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4. G24.78+0.08 

G24.78+0.08 is a high-mass star forming region that con- 
tains a cluster of massive YSOs in different evolutionary 
stages. Furuya et al. (|2002[1 firstly detected four molecu- 
lar cores in the region, identified as A, B, C, and D in 
Fig. HI U| l. with two of them, core A and C, powering a 
bipolar molecular outflow. When observed with higher an- 
gular resolution, the continuum and methyl cyanide emis- 
sion toward core A is resolved into two separate cores, 
named Al and A2 in Paper I. It should be mentioned 
that although one cannot rule out the possibility that the 
continuum and methyl cyanide peaks in A2 could be the 
result of the interaction of the outflow with the surround- 
ing dense gas, we consider them as two distinct cores in 
our study (see a more detailed discussion in Sects. 14.2.21 
and l4~!T3jl . 

4.1. Results 

In this section we analyze the line and continuum observa- 
tions carried out with the BIMA and PdBI interferometers 
toward G24.78+0.08. Subsequently to the lower angular 
resolution BIMA data, we present the higher angular res- 
olution PdBI data. By presenting the results in such a 
way, the analysis of G24.78+0.08 starts with the study of 
the large scale molecular clumps, and then continues with 
the study of the smaller scale embedded cores. 

4.1.1. N 2 H+, C 34 S and CH 3 OH 

N 2 H+ (1-0), C 34 S (2-1), and CH 3 OH (2_i-l_i) E, 
(2o-lo) A+, (2o-lo) E were observed toward the hot core 
G24.78+0.08 with the BIMA interferometer. Figure EH 
shows the integrated emission around the systemic veloc- 
ity of 111 km s" 1 , overlaid on the 3.16 mm continuum 
emission, observed also with BIMA. As may be seen in 
this figure, a well defined molecular core is found in asso- 
ciation with G24.78+0.08, although the low angular res- 
olution of the BIMA maps (8 ~ 14") makes it impossi- 
ble to disentangle the emission of each core identified in 
the region by Furuya et al. I|2002|l . This molecular clump 
has also been detected in several rotational transitions of 
CO and its isotopomers, and in CS (3-2) (Cesaroni et 
al. |5003 ) . In Fig. EH one can see that each species peaks 
at different positions: C 34 S toward core A, and CH3OH 
and N 2 H + between cores C and A, with the N2H + emis- 
sion peaking much closer to core C. This suggests that ev- 
ery species could be tracing different physical conditions 
and/or chemical evolution of the molecular gas toward the 
YSOs. The lower resolution of our observations, however, 
did not allow us to investigate this in more detail. Taking 
into account that we have detected rotating structures as- 
sociated with some of the individual cores embedded in 
G24.78+0.08 (see Paper I), we also searched for the pres- 
ence of a global velocity gradient that could be indicative 
of large scale rotation of the whole cloud. To do this we 
fitted the C 34 S (2-1) lines with a Gaussian at each point 



G24.78 + 0.08 - N 2 H + (1-0) - 12 CO(1-0) 




a(J2000) 

Fig. 11. Overlay of the BIMA N 2 H+ (1-0) integrated 
emission (greyscale), and the 12 CO (1-0) integrated emis- 
sion obtained with the IRAM 30-m telescope by Cesaroni 
et al. lj271T)3l contours) toward G24.78+0.08. The 12 CO 
emission has been integrated under the line, and the 
contour levels range from 4 to 194 Kkms^ 1 in steps 
of 10 Kkm s _1 . Contours for N 2 H + are the same as in 
Fig. The triangle pinpoints the position of the peak of 
the 20 cm continuum emission extracted from the NVSS 
database. The dashed circumference represents the BIMA 
primary beam (50% attenuation level). 



in order to measure shifts from the systemic velocity of 
the cloud, and we also minimized the \ 2 measure of the 
difference between a model with a linear velocity gradient 
and the C 34 S (2-1) data at each point. However, no sig- 
nificant rotation was detected by any of the two methods, 
except for the central part of the cloud where the velocity 
gradient is more evident and consistent with that found 
on a smaller scale in CH3CN. 

A second molecular clump has clearly been detected 
in N 2 H + and is also visible in CH3OH and marginally in 
C 34 S, toward the northwest of G24.78+0.08 (see Fig.[TO|). 
Millimeter continuum emission has also been detected in 
our 3.16 mm low angular resolution (9 ~ 18") BIMA ob- 
servations toward the northwest. This continuum emission 
peaks at the position of an extended H11 region visible at 
20 cm in the NVSS survey (Condon et al. 11998(1 . and is 
probably associated with it. However, as can be seen in 
the figure, there is also some continuum emission associ- 
ated with the N 2 H + clump. Submillimeter emission asso- 
ciated with this clump has also been mapped at 850 and 
450/im with the Submillimeter Common-User Bolometric 
Array (SCUBA) by Walsh et al. PTDljl . In their maps 
it is also visible lower-level extended emission connecting 
the two clumps in the region. The clump has also been 
detected in several rotational transitions of CO and its 
isotopomers (Cesaroni et al. 2003). However, the peak of 
the 12 CO emission is displaced from the position of the 
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Fig. 12. Combined BIMA+PdBI map of the 3.3 mm continuum emission (left panel), and PdBI map of the 1.4 mm 
continuum emission (right panel) continuum emission toward the core of G24.78+0.08. The contour levels range from 
2 to 20 mJy beam -1 in steps of 2 mJy beam -1 and from 20 to 80 mJy beam -1 in steps of 10 mJy beam -1 for the 
3.3 mm map, and from 20 to 200 mJy beam -1 in steps of 60 mJy beam -1 for the 1.4 mm map. The synthesized beam 
is shown in the lower left-hand corner. The crosses mark the position of the YSOs identified in the region (Furuya et 
al.EEEl Paper I; see Sect. 11X2) . 



N2H + emission peak (see Fig, lilt , unlike the molecular 
clump at the center of the map for which both CO and 
N2H + peak at the same position. The lack of CO emission 
toward this northern clump could be due to an opacity ef- 
fect or to molecular depletion. For low-mass starless cores, 
the abundance of some molecules, including CO and CS, 
dramatically decreases toward the core center, with re- 
spect to the outer region. On the other hand, molecules 
such as N2H + are found to have constant abundances (e.g. 
Tafalla et al. 2002 ) . These results are explained by the de- 
pletion of carbon and sulfur bearing molecules onto dust 
grains at high densities (a few times 10 4 cm -3 ) and low 
temperatures occurring in dense core interiors (e.g. Bergin 
& Laneer Il997|) . Nitrogen bearing molecules are instead 
unaffected by this process up to densities of several times 
10 5 cm -3 , or even 10 6 cm -3 . In order to assess that de- 
pletion is the cause of the lack of CO emission we need to 
estimate the opacity and the temperature of this northern 
clump. However, only the CO (1-0) and (2-1) lines have 
been observed toward this clump: therefore additional ob- 
servations of an optically thin CO isotopomer and of a 
good temperature tracer are needed. Unfortunately, the 
energy of the three CH3OH (2^-1^) transitions observed 
with BIMA are too similar to derive any temperature by 
means of the RD method. In case that the CO freeze-out 
is confirmed, the northern clump would represent the first 
example of a chemically depleted high-mass clump. 

4.1.2. Continuum emission 

In Fig. El we show the maps of the combined BIMA + 
PdBI 3.3 mm continuum emission and the PdBI 1.4 mm 
continuum emission toward G24.78+0.08. As mentioned 
in Sect. |4~1~T1 3.16 mm continuum low angular resolution 



(6 ~ 18") BIMA maps are shown in Fig.HHI The highest 
angular resolution (9 ~ 1") observations at 3.3 mm have 
resolved the continuum emission into three sources, A, 
B, and C, already detected at millimeter wavelengths by 
Furuya et al. I)2002|l . The continuum emission at 3.3 mm is 
dominated by an extended structure surrounding cores A 
and B. The strongest core detected at 3.3 mm wavelength 
in the region is core A. The emission clearly shows two 
components, a centrally peaked source, plus an extended 
component elongated toward the northwest. In fact, the 
1.4 mm dust emission associated with core A has been 
resolved into two separate cores, cores Al and A2 (see 
Fig. 1121 and Paper I), with core A2 located in the direc- 
tion of the elongation toward NW seen at 3.3 mm. The 
two cores are separated by ~ l'/5 or ~ 11500 AU. The 
positions, fluxes and deconvolved diameters of the cores 
are given in Table |21 The other core detected at both 
wavelengths is C. The continuum emission associated with 
this YSO is much less extended, and it has only been re- 
solved at 1.4 mm. The measured flux densities of cores A 
and C at 3.3 mm are consistent with the previous mea- 
surements by Furuya et al. (|2()f)2l see their Fig. 3). At 
1.4 mm, the sum of the flux densities of cores Al and 
A2 is consistent with the expected flux density for core 
A, while for C the flux that we have measured is slightly 
lower, but still consistent, with the model predictions of 
Cesaroni et al. (2003). The third source detected in the 
region at 3.3 mm is the UC Hn region B (Codella et 
al. 119971 see Fig Wl\ , which is slightly elongated in the 
southeast-northwest direction. This has not been detected 
at 1.4 mm; however, the 1.4 mm flux expected on the basis 
of the results of Cesaroni et al. ( 2003 ) is lower than the 3cr 
noise level of our map, which is about 20 mJy beam -1 . 
The flux at 3.3 mm is also consistent with the values of 
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Fig. 13. Same as Fig.0for cores A1+A2 in G24.78+0.08. 



Cesaroni et al. ( 2003). Finally, we have not detected at 
any wavelength core D, which has been detected at 2.6 
and 2 mm by Furuya et al. (2002). However, the expected 
flux of this core at 3.3 and 1.4 mm, according to the pre- 
dictions of Cesaroni et al. fZOlQ 3j , falls below the 3cr noise 
level of our maps, 2 and 20 mJy, respectively. 

4.1.3. CH 3 CN and CH 3 13 CN 

PdBI spectra of the CH 3 CN (12-11) and (5-4), CH 3 13 CN 
(12-11) and (5-4), and CH 3 CN (5-4) v s = 1 lines toward 
cores Al and A2 in G24.78+0.08 are shown in Fig.lTSl The 
spectra have been obtained averaging the emission over 
the 3cr contour level area, which surrounds both cores Al 
and A2. Figure HI shows the PdBI maps of the CH 3 CN 
(12-11) emission averaged under the K =1 and 2 compo- 
nents, and under the K = 6, 7 and 8 components, and the 
CH 3 13 CN emission averaged under the K =1 and 2 com- 
ponents. Figure [TBI shows the maps of the CH 3 CN (5-4) 
emission averaged under the K — to 4 components, the 
CH 3 13 CN emission averaged under the K = to 3 com- 
ponents, and the CH 3 CN vg = 1 emission averaged under 
all observed transitions. 

The methyl cyanide emission of the two cores has not 
been resolved at 3.3 mm because the separation between 
them is smaller than the synthesized beam, whereas they 
have been resolved at 1.4 mm (Figs. ITHand ll5|l . For the 
sake of comparison with the CH 3 CN (5-4) lines, also for 
the CH 3 CN (12-11) emission we decided to plot the spec- 
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Fig. 14. PdBI intensity maps of the CH 3 CN (12-11) emis- 
sion averaged under the K = 1 and 2 components (top 
panel), under the K = 6, 7, and 8 components [middle 
panel), and of the CH 3 13 CN (12-11) emission under the 
K = 1 and 2 components (bottom panel) toward cores Al 
and A2 in G24.78+0.08. The contour levels range from 
0.14 to 0.98 Jy beam -1 (top panel), and to 0.56 Jy beam -1 
(middle and bottom panels), in steps of 0.14 Jy beam -1 . 
The synthesized beam is shown in the lower left-hand cor- 
ner. The white crosses denote the position of the 1.4 mm 
continuum emission peaks. The triangle marks the posi- 
tion of the VLA 1.3 cm continuum emission peak (Codella 
et al. lTM7| . 

tra integrated over a region surrounding both cores. As 
can be seen in Fig. 1131 several if-components of the differ- 
ent rotational transitions of methyl cyanide were clearly 
detected. No methyl cyanide emission has been detected 
toward cores B or D with the PdBI. For core C, only the 
CH 3 CN (5-4) emission has been barely detected, with an 
averaged intensity peak of ~ 30 mJy beam -1 . Cores B 
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Fig. 15. PdBI intensity maps of the CH 3 CN (5-4) (top 
panel) emission averaged under the K = to 4 com- 
ponents, the CH3 13 CN (middle panel) emission averaged 
under the K = to 3 components, and of the v$ = 1 emis- 
sion (bottom panel) averaged under the (1,1), (2,1), (0,1), 
(3,1), (1,-1), (4,1), (2,-1), (3,-1), (4,-1), and (1,1) tran- 
sitions toward cores Al and A2 in G24.78+0.08. The con- 
tour levels range from 0.05 to 0.35 Jy beam -1 in steps of 
0.05 Jy beam -1 (top panel), from 0.01 to 0.06 Jy beam -1 
in steps of 0.01 Jy beam -1 (middle panel), and from 0.005 
to 0.030 Jy beam -1 in steps of 0.005 Jy beam -1 (bottom 
panel). The synthesized beam is shown in the lower left- 
hand corner. The white crosses denote the position of the 
1.4 mm continuum emission peaks. The triangle marks 
the position of the VLA 1.3 cm continuum emission peak 
(Codella et al. lTOI . 



and C have been detected in the CH3CN (8-7) transition 
by Furuya et al. Ij2002|l . hence the fact that they are not 
detected in the (12-11) lines could be due to these cores 
being significantly colder than cores Al and A2. As done 
for the spectra toward G31.41+0.31, we fitted the multi- 
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Fig. 16. BIMA integrated intensity map of the CH 3 CN 
(5-4) emission under the K = 0, 1, 2, and 3 components 
toward cores Al and A2 in G24.78+0.08. The contour lev- 
els are 3, 6, 9, 12, 15, and 18 times 0.35 Jy beam -1 km s -1 . 
The synthesized beam is shown in the lower left-hand cor- 
ner. The white crosses mark the position of the 1.4 mm 
continuum emission peaks. 



pie if-components of the differ ent CH 3 CN and CH 3 13 CN 
rotational transitions (see Sect. l3.1?2l for details on the fit- 
ting). The CH 3 CN and CH 3 13 CN (12-11) line parameters 
for cores Al and A2, individually, are given in Tableland 
El respectively, while the CH 3 CN and CH 3 13 CN (5-4), and 
CH 3 CN (5-4) Ug = 1 line parameters for the emission sur- 
rounding both cores Al and A2 (hereafter A1+A2) are 
given in Table l9"l and ITOI As can be seen in Tables Hand 
03 the systemic velocities and the line widths derived for 
the (12-11) transition are quite similar for both cores: 
Vlsr ~ 110.6-110.9 km s -1 , and FWHM ~ 7 km s -1 . 
The systemic velocity is also consistent with that found 
for the ground state and v$ = 1 (5-4) transitions in cores 
A1+A2, although those derived from the ground state 
lines are slightly higher (111.2-111.6 km s -1 ). Regarding 
the line widths, the values derived for the (5-4) transitions 
are quite different, with ~ 7.8 km s -1 for CH 3 CN and 
~ 4.6 km s -1 for CH 3 13 CN. As done for G31.41+0.31 (see 
Sect. 13. Q|| . one can derive the ratio between the observed 
and the intrinsic line widths, which for G24.78+0.08 is 
1.7 ± 0.04. On the other hand, the ratio derived from the 
optical depths of the CH 3 CN and CH 3 13 CN (5-4) lines 
(assuming a relative abundance [CH 3 CN]/[CH 3 13 CN]=40; 



see Wilson & Rood 1994) ranges from 1.5 to 1.9. Thus, also 
in this case, optical depth effects could be responsible for 
the variation of the line widths. The line widths found for 
(5-4) ug = 1 are slightly smaller but still consistent with 
those found for the (12-11) transition. 

Both core Al and A2 are clearly visible in Fig. El be- 
ing the emission associated with core A2 the strongest 
for all the transitions. The emission of both cores has 
not been resolved at 3 mm (Fig. I15|) . As can be seen in 
Fig. 1141 the morphology of the emission is quite similar 
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Fig. 17. Overlay of the integrated intensity map of the 
CH 3 CN (5-4) emission under the K = 0, 1, and 2 
components obtained with the PdBI {contours) on the 
map obtained with BIMA (greyscale). The contour lev- 
els range from 0.05 to 0.35 Jy beam -1 km s _1 (greyscale) 
and to 0.40 Jy beam _1 km s _1 (contours) in steps of 
0.05 Jy beam" 1 km s" 1 . The BIMA map has been ob- 
tained by restoring the 'clean' map with the same beam 
than that of the PdBI observations (see text). The syn- 
thesized beam is shown in the lower left-hand corner. The 
white crosses mark the position of the 1.4 mm continuum 
emission peaks. 

G24. 78 + 0. 08 CORES A1&A2 BIMA & PdBI CH 3 CN(5-4) 




Fig. 18. Comparison between the BIMA (thick line) and 
the PdBI (thin line) CH3CN (5-4) spectra toward cores 
A1+A2 in G24.78+0.08. Both spectra have been obtained 
averaging the emission over the same region. The BIMA 
channel maps have been obtained by restoring the 'clean' 
map with the same beam than that of the PdBI ob- 
servations (see text). The number of the different K- 
components are marked with dashed lines in the upper 
part. 



for the different if-components of CH 3 CN and CH 3 13 CN: 
one sees two compact cores, roughly coincident with the 
two 1.4 mm continuum cores Al and A2, plus an extended 
component. Such a good agreement between the gas and 
the dust emission indicates that they are tracing the same 
material. In addition, there is no clear evidence of any 



emission dip for lower isT-components of CH3CN toward 
the position of the embedded sources in the two cores, as 
for G31.41+0.31. This suggests that the temperature gra- 
dient in cores Al and A2 in G24. 784-0. 08 is less prominent 
than in G31.4l4-0.31. This is confirmed by the T, ■ t mea- 
sured in each core, as illustrated in Sect. I4.2T1 In Fig. 1141 
we also mark the position of the peak of the 1.3 cm contin- 
uum emission of the UC Hn region detected with the VLA 
by Codella et al. {1997). Note that this H11 region is as- 
sociated with core Al, being located at ~ O'.'l (~ 800 AU 
at the distance of the region) of the millimeter continuum 
source, whereas no free-free emission is reported toward 
core A2. 

Figure [H shows the BIMA integrated CH 3 CN (5-4) 
emission under the K = 0, 1, 2, and 3 components toward 
cores Al and A2. The cores have been clearly resolved in 
these observations (0 ~ Of! 7). As can be seen in the fig- 
ure the strongest emission comes from core A2, with the 
peak slightly displaced from that of the 1.4 mm continuum 
emission. On the other hand, the emission from core Al is 
hardly detected. This could be due to a temperature effect, 
being core Al probably hotter than core A2. We checked 
whether it could be an opacity effect by computing the 
optical depth in the cores from the ratio of CH 3 13 CN 
and CH 3 CN (12-11) K = 2 component (see Sect. EEEfl . 
However, the optical depths obtained are similar for both 
cores, Ti2 ~ 12, suggesting that opacity cannot account 
for the difference in CH 3 CN (5-4) emission between both 
cores. In Fig. El we compare the CH 3 CN (5-4) maps ob- 
tained with BIMA and the PdBI. To compare them we 
have first reduced the resolution of the BIMA map by 
restoring the 'clean' BIMA map with the same beam as 
the PdBI observations. As can be seen in this figure there 
is an excellent agreement between the two maps, which 
confirms that most of the CH 3 CN (5-4) emission is asso- 
ciated with core A2. Note that some extended emission is 
not shown by the BIMA map, probably due to the fact 
that it is filtered out by the interferometer due to the lack 
of short spacings. In Fig.^lthe same comparison is shown 
for the spectra obtained integrating the emission over the 
same region. Again, taking into account that some of the 
emission is probably filtered out by the BIMA interferom- 
eter, one can say that the agreement between both spectra 
is remarkable. 



4.2. Discussion 

4.2.1. Temperature and column density 

We derived the rotational temperature and total col- 
umn density of the methyl cyanide molecules by means 
of the RD method (see Sect. I3.1.2|) for cores A14-A2 in 
G24.784-0.08 (see Fig. "TJJ. The insufficient angular res- 
olution of the 3 mm observations makes unable to sepa- 
rate the emission from the two cores, and thus the emis- 
sion has been integrated over an area of ~ 4" of diam- 
eter (~ 12 arcsec 2 ) that covers both core Al and A2. 
As stated before, such a region would correspond to core 
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Table 7. CH 3 CN and CH 3 13 CN (12-11) line parameters for core Al in G24.78+0.08 



Line 


Vlsr 


FWHM 


n 


J'T B dV 


rms 


K 


(km s" 1 ) 


(km s^ 1 ) 


(K) 


(K km s" 1 ) 


(K km s" 1 ) 


CH 3 CN 





110.82±0.02 


7.25 ± 0.03 


26.8 ± 1.5 


207 ±3 


0.9 


1 


)) 


tt 


27.6 ±1.5 


213 ±3 


tt 


2 


tt 


tt 


28.3 ± 1.5 


219 ± 3 


tt 


3 


tt 


tt 


30.3 ±1.5 


234 ±2 


tt 


4 


110.82±0.05 


7.12 ±0.07 


23.6 ± 1.3 


179 ±3 


1.0 


5 




tt 


18.2 ± 1.3 


207 ±3 


55 


6 


110.64 ± 0.12 


6.90 ±0.16 


17.0 ±0.8 


125 ± 5 


0.8 


7 


J) 




7.7 ±0.8 


56 ±3 


0.8 


8 


tt 


tt 


3.4 ±0.8 


25 ±3 


tt 


9 


?) 




(b) 


(b) 


tt 


10 


tt 


55 


< 0.8 




tt 


11 


)5 




< 0.8 




tt 


CH 3 13 CN 





110.82±0.05 


7.12 ±0.07 


(c) 


(c) 


1.0 


1 


?! 


5? 


7.0 ± 1.3 


53 ±3 


tt 


2 


J) 


tt 


6.6 ± 1.3 


50 ±2 


tt 


3 


110.64 ±0.12 


6.90 ±0.16 


8.7 ±0.8 


64 ±5 


0.8 


4 


110.64 ±0.10 


6.92 ±0.14 


2.6 ±0.8 


19 ±3 


0.7 


5 


tt 


It 


0.9±0.8 d 


6±3 d 


tt 


6 


110.64 ±0.12 


6.90 ±0.16 


1.8 ±0.8 


13 ±3 


0.8 


7 


110.68 ±0.72 


7.0 C 


< 0.8 




0.8 


8 


tt 


tt 


3.4±0.5 f 


26±5 f 


tt 


9 


?) 


5? 


1.3 ±0.5 


10 ±5 


tt 



(a) Brightness temperature integrated over a 3cr emission level area around core Al. 

(b) Blended with 13 CO (2-1). 

(c) Blended with CH 3 CN K = 5. 

(d) Blended with CH 3 CN K = 7. 

(e) This parameter was held fixed in the fit. 

(f) Possibly blended with wing emission of 13 CO (2-1). 
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Fig. 19. 

G24.784 



Same 
0.08. 



as Fig. El for the cores A1±A2 in 



A of Codella et al. fiWfy and Furuya et al. ( 2002). As 
already done for G31.41+0.31, we used only the ground 
state transitions of the isotopomer CH3 13 CN, and those 
from the CH3CN v$ = 1 state in the fit performed to the 
Boltzmann plot, because the ground level transitions (see 
Fig. I13|) appear to be optically thick. We assumed a rel- 
ative abund ance [CH 3 CN]/[CH 3 13 CN]=40 (see Wilson & 
Rood ■ The T rot and source averaged iVcH 3 CN that 
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Fig. 20. Rotation diagram of CH 3 13 CN (12-11) for cores 



Al (top) and A2 (bottom) in G24.78- 
posed fits. 



0.08, with supcrim- 
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Table 8. CH 3 CN and CH 3 13 CN (12-11) line parameters for core A2 in G24.78+0.08 



Line 


Vlsr 


FWHM 


n 


J'T B dV 


rms 


K 


(km s" 1 ) 


(km s" 1 ) 


(K) 


(K km s" 1 ) 


(K km s" 1 ) 


CH 3 CN 





110.68±0.03 


7.25 ± 0.03 


22.0 ± 1.6 


170 ±3 


1.1 


1 


11 


n 


22.6 ± 1.6 


174 ±3 


11 


2 


11 


11 


23.8 ± 1.6 


184 ±3 


11 


3 


11 


11 


26.2 ± 1.6 


202 ± 2 


11 


4 


110.59±0.07 


7.18 ±0.10 


20.8 ±0.9 


159 ±3 


1.3 


5 






16.0 ±0.9 


123 ± 5 


11 


6 


110.87 ±0.08 


6.74 ±0.11 


16.3 ± 1.1 


117 ±3 


0.6 


7 


11 


J) 


7.9 ± 1.1 


57 ±2 


11 


8 


110.88 ±0.10 


6.66 ±0.13 


4.3 ± 1.0 


30 ±3 


0.7 


9 


110.50 ±0.60 


7.0 b 


(c) 


(c) 


0.8 


10 


11 


)} 


< 0.8 




)) 


11 


11 




< 0.8 




11 


CH 3 i3 CN 





110.59±0.07 


7.18 ±0.10 


(d) 


(d) 


1.3 


1 


11 




6.1 ±0.9 


47 ±3 


11 


2 


?) 




6.0 ±0.9 


46 ±3 


11 


3 


110.87 ±0.08 


6.74 ±0.11 


9.3 ± 1.1 


67 ±3 


0.6 


4 


?) 


n 


3.3 ± 1.1 


24 ±2 


11 


5 


11 


11 


1.8 ± 1.1° 


13±3 C 


11 


6 


110.88 ±0.10 


6.66 ±0.13 


2.1 ± 1.0 


15 ±3 


0.7 


7 


110.50 ±0.60 


7.0 b 


< 0.8 




0.8 


8 


11 


)} 


3.7±0.6 f 


27±5 f 


11 


9 


5) 




1.2 ±0.6 


9±5 


11 



(a) Brightness temperature integrated over a 3u emission level area around core A2. 

(b) This parameter was held fixed in the fit. 

(c) Blended with 13 CO (2-1). 

(d) Blended with CH 3 CN K = 5. 

(e) Blended with CH 3 CN K = 7. 

(f) Dominated by wing emission of 13 CO (2-1). 

Table 9. CH 3 CN and CH 3 13 CN (5-4) line parameters for core Al±A2 in G24.78±0.08' 



Line 


Vlsr 


FWHM 




J'T B dV 


rms 


K 


(km s- 1 ) 


(km s- 1 ) 


(K) 


(K km s- 1 ) 


(K km s- 1 ) 


CH 3 CN 





111.61 ±0.02 


7.74 ± 0.03 


10.6 ± 1.0 


87 ±1 


0.3 


1 


11 




8.9 ± 1.0 


70 ±1 


)) 


2 


11 


11 


10.8 ± 1.0 


89 ±1 


;j 


3 


11 


11 


11.6 ± 1.0 


96 ±1 


)) 


4 


11 


11 


4.9 ± 1.0 


40 ± 1 


)? 


CH 3 "CN 





111.20 ±0.08 


4.58 ±0.09 


2.00 ± 0.27 


9.8 ±0.4 


0.18 


1 


)) 


11 


1.51 ±0.27 


7.39 ±0.35 


)) 


2 


j; 


11 


0.86 ±0.27 


4.21 ±0.35 




3 


)? 


11 


1.33 ±0.27 


6.50 ± 0.35 


)? 


(a) Emission integrated 


over a region 


surrounding 


cores Al and A2 


in G24.78±0. 



(b) Brightness temperature integrated over a 3<r emission level area. 



we derived by taking into account all transitions are 132 K 
and 2 x 10 16 cm~ 2 , respectively. The T ro t and A^cHaCN 
that we estimated from each set of J ± 1 — > J transitions 
of CH 3 13 CN are 102 ± 5 and 130 ± 2 K, and 1 x 10 16 and 
2 x 10 16 cm~ 2 for the (5-4) and (12-11) transitions, re- 
spectively. Note that unlike the case of G31.41±0.31, for 
G24.78±0.08 these two estimates are quite similar, which 



could be indicating that all the transitions are tracing 
more or less the same gas or, as mentioned before, that the 
temperature gradient is less evident. From NH 3 emission 
averaged over a similar area Codella et al. 1|1997|) have de- 
rived a slightly lower kinetic temperature, Tk = 87 K: this 
is not surprising because NH 3 is probably tracing a differ- 
ent region, somewhat more extended and colder. The es- 
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Table 10. CH 3 CN (5-4) v 8 = 1 line parameters for core A1+A2 in G24.78+0.08 a 



Line 


Vlsr 
(km s \ 


1 


FWHM 

fi 1 \ 

(km s 1 ) 


(K) 


j T B dV 
(K km s" 1 ) 


rms 
(K km s _1 ) 


(1,1) 


110.61 ±0. 


18 


6.26 ±0.19 


0.77 ±0.14 


5.15 ±0.34 


0.07 


(2,1) 


110.84 ±0. 


12 


5.78 ±0.13 


0.82 ±0.14 


5.04 ±0.32 


0.07 


(0,1) 


)) 




it 


0.91 ±0.14 


5.61 ±0.32 




(3,1) 


)) 






0.22 ±0.14 


1.38 ±0.31 


)} 


(1,-1) 


)) 






0.83 ±0.14 


5.11 ± 0.32 


)} 


(4,1) 


110.61 ±0. 


18 


6.26 ±0.19 


0.60 ±0.14 


3.98 ±0.34 


0.07 


(2,-1) 


)) 






0.73 ±0.14 


4.84 ± 0.34 




(3,-1) 


V 






< 0.07 






(4,-1) 


110.61 ±0. 


18 


6.26 ±0.19 


0.11 ±0.14 


0.71 ±0.32 


0.07 


(1,1) 


110.61 c 




6.26 c 


0.89 ± 0.09 d 


5.94±0.52 c 


0.04 



(a) Emission integrated over a region surrounding G24.78±0.08 cores Al and A2. 

(b) Brightness temperature integrated over a 3<r emission level area. 

(c) This parameter was held fixed in the fit. 

(d) Computed from the integrated line intensity. 

(e) Measured with the command print area of CLASS. 



timate of the gas column density derived by these authors 
is consistent with ours. Taking into account that with the 
angular resolution achieved at 1 mm it is possible to sepa- 
rate the emission from cores Al and A2, we made rotation 
diagrams for the CH3 13 CN (12-11) for each of them sep- 
arately (see Fig. I20|) . As seen in this plot, the T rot and 
source averaged -/Vch 3 cn are similar for both cores. 

In Fig. we show the T rot map and the Nqh 3 cn map 
overlaid to the CH3 13 CN (12-11) emission averaged un- 
der the K =1 and 2 components toward cores Al and 
A2. We computed the maps as done for G31.41±0.31 (see 
Sect. 13.1.2(1 . This figure shows that both cores have sim- 
ilar T rot , ~ 120-150 K, consistent with the source aver- 
aged temperatures derived in the previous paragraph. In 
Sect. l4.1^1 we proposed that core Al could be hotter than 
A2, because the former is hardly detected in CH3CN (5-4) 
(see Figs.EDandEIl. The fact that we derive similar tem- 
peratures for both cores by means of the RD method rules 
out this hypothesis and indicates that the assumption of 
LTE is not valid in this case. In order to explain the dif- 
ference in CH3CN (5-4) emission, one should take into 
account the IR radiation field of the UC H11 region asso- 
ciated with this core (see Fig. 114(1 . Figure PHI also shows 
that there is no clear radial temperature variation in the 
cores, unlike for G31.41±0.31 (see Fig. 01. Regarding the 
column density map, -/Vch 3 cn, one can see that there are 
two maxima roughly at the position of the methyl cyanide 
and millimeter continuum emission peaks for cores Al and 
A2 (Fig l21(l . confirming that for these two cores, CH3CN 
and dust trace the same material. 

4.2.2. Nature of core A and evolutionary sequence 

Core A was first detected at 1.3 cm by Codella et al. ( 1997) 
as a very compact and dense unresolved NH3 core. The 
gas and dust millimeter emission from this core is the 
strongest and most extended of all the cores detected in 
the region. In addition, Furuya et al. 12)02) have mapped 



a bipolar molecular outflow that is probably powered by a 
YSO embedded in the core. As already mentioned, thanks 
to the higher angular resolution of our PdBI millimeter 
observations, we have resolved core A into two separate 
cores, Al and A2 (Paper I). These have been detected in 
gas and dust emission and both could be hosting high- 
mass YSOs. In such a case, the YSOs might be in dif- 
ferent evolutionary stages, because the two cores hosting 
them have roughly the same mass (~ 130 M© for Al and 
~ 80 Mq for A2; see Paper I), but only Al has been found 
in association with an UC Hn region, i.e. a more evolved 
YSO (see Fig.IUJ. The UC Hn region has a diameter (at 
50% of the peak emission) of ~ 0715 or — 1200 AU and is 
hence very small and probably in an early stage of its evo- 
lution. Regarding the bipolar outflow detected by Furuya 
et al. 1(2002(1 in association with A, it is not evident which 
of the cores, Al or A2, hosts the powering source. Core Al 
lies slightly closer to the geometrical center of the outflow, 
but the small separation between the cores and the fact 
that they are aligned along the outflow axis make it diffi- 
cult to establish any association unambiguously. Sensitive 
mid-IR observations would be very helpful to assess the 
YSO content of each core and possibly constrain their evo- 
lutionary state. 



Furuya et al. I2 0f)2^l have concluded that the hot core 
G24.78±0.08 contains at least 4 high-mass YSOs with dif- 
ferent ages: more precisely, £b > > tc > tu- I n the 
light of our results, we can now infer that the embed- 
ded YSOs could be 5, being Al probably older than A2. 
In addition, the fact that source C is colder (T ~ 30 K; 
see Codella et al. HW7\ than source A2 (T ~ 127 K; see 
Sect. 14.2.1(1 suggests that A2 is probably older than C. 
Thus, the ages of the YSOs in G24.78±0.08 seem to be in 
the order <b > £ai > ^A2 > tc > tn- 
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Fig. 21. Overlay of the PdBI CH 3 13 CN (12-11) emission 
averaged under the K = 1 and 2 components {contours) 
on the T rot map (top panel) and the -/Vch 3 cn map {bottom 
panel) in greyscale, derived by fitting the Boltzmann plots 
at each position, toward cores Al and A2 in G24.78+0.08. 
The contour levels range from 0.14 to 0.56 Jy beam -1 in 
steps of 0.56 Jy beam -1 . Greyscale levels range from 60 
to 180 K by 30 K {top panel), and from 7 x 10 15 to by 
1.7x 10 17 cm" 2 times by 2x 10 16 cm -2 {bottom panel). The 
white crosses mark the position of the 1.4 mm continuum 
emission peaks. 



4.2.3. Velocity field: the rotating toroid 

In Paper I we report on the detection of a clear veloc- 
ity gradient in both core Al and A2, with Vlsr increasing 
steadily along a direction at a P. A. ~ 50°, which is perpen- 
dicular to the axis of the molecular outflow detected to- 
ward core A by Furuya et al. (jj002 ). As already discussed 
in Paper I, the most plausible interpretation for such ve- 
locity gradients is that the cores have toroidal structures 
undergoing rotation about the corresponding outflow axis. 
In such a scenario, the two cores Al and A2 would corre- 
spond to two distinct rotating toroids. However, the sep- 
aration of the cores is small, l'/5 or 12000 AU, and we 
cannot rule out the possibility that they are part of the 
same geometrically thick rotating toroid, or pipe, elon- 
gated along the direction of the outflow. In such a sce- 
nario, assuming core Al as the powering source of the 
outflow, the velocity gradient detected in the region could 
be tracing a general rotation of this pipe, related to the 
rotation associated with the launching of the jet. Taking 
into account jet models that predict that the flow moves 
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Fig. 22. Same as Fig.[H]along the major axis of the toroid 
(P.A. - 50°) in core Al in G24.78+0.08. The systemic 
velocity of G24.78+0.08 is 111 km s -1 . The contour levels 
are 4, 12, 20, 28, 36, 44, and 50 Jy beam" 1 . 



along the magnetic surfaces conserving its angular mo- 
mentum after the Alfven surface (e.g. Spruit IT996I Konigl 
& Pudritz 2000), the jet launching radius should be about 
half of the radius of the toroid in Al in order to reach 
the maximum velocity in A2, which is measured at a dis- 
tance of ~ 0.02 pc. Such a launching radius is consistent 
with the data, and thus, it calls for further higher angular 
resolution observations aimed to face the open question. 

Due to the difficulty to establish which of the two cores, 
Al or A2, is powering the bipolar outflow, we arbitrarily 
assume that it is core Al, as also done in Paper I, be- 
cause its association with an UC Hn region indicates that 
a YSO has already formed. Thus, we model only the rotat- 
ing structure associated with core Al. Note that the PV 
plots for both cores are very similar, and thus, the results 
and conclusions derived from the modeling, with the ex- 
ception of the exact value of the model parameters, apply 
also to core A2. The velocity shift measured over an extent 
of ~ 9000 AU is ~ 3 km s _1 , that is, a velocity gradient of 
about 70 km s _1 pc -1 . The kinematics of the gas toward 
core Al can be investigated with a PV cut of the CH3CN 
(12-11) K = 3 emission along a direction (P.A. ~ 50°) 
roughly perpendicular to the outflow axis (Fig. I22f> . As 
done for G31.41+0.31, we modeled the emission by assum- 
ing power-law distributions for the temperature, and the 
velocity field in the toroid (see Sect. I3.2.2JI . and adopted 
constant density in the core. The models computed were 
with constant rotation velocity and with constant angu- 
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Fig. 23. Line width plot of the CH 3 CN (12-11) toward 
core Al in G24.78+0.08. The contour levels range from 5.5 
to 7.5 km s -1 in steps of 0.5 km s _1 . The crosses mark the 
positions of the 1.4 mm continuum emission peaks. The 
white triangle marks the position of the 1.3 cm continuum 
emission peak from the Hll region (Codella et al. 11997) . 



lar velocity, because Keplerian rotation is not possible 
due to the huge mass of the toroid. For the temperature 
distribution we again adopted a temperature power law 
T oc i?~ 3 / 4 . As can be seen in Fig. 1221 from the modeling 
itself it is not possible to distinguish between constant an- 
gular velocity (v oc i?; top panel) or constant rotation (v = 
constant; bottom panel), because both synthetic PV dia- 
grams fit the data reasonably well. We computed models 
with and without temperature gradients, and both fitted 
the data equally well. The fact that a molecular outflow 
and an UC Hll region have been detected in association 
with this core implies the existence of an embedded source 
powering it. Therefore, although not clearly seen in our 
data, one may reasonably assume the existence of a tem- 
perature gradient driven by the central heating source. 
Thus, we adopted a temperature distribution T oc i?~ 3 / 4 
for the models. The fact that the G24.78+0.08 CH 3 CN 
emission, unlike that of G31.41+0.31, has not got a "hole" 
at the center, is consistent with the gas emission of core 
Al in G24.78+0.08 being optically thinner than that of 
G31.41+0.31. 

The best fit model parameters are given in Table H3 
The value of R inn is 0"3 (~ 2300 AU), which is the size 
of the UC Hll region at 1.3 cm at the 3cr emission level 
as recently observed with the VLA in the A configura- 
tion (Beltran, private communication). The value of v TO t 
is consistent with that derived directly from the velocity 
gradient in Paper I. The i? ut derived from the model is 
roughly twice the value measured from the 1.4 mm contin- 
uum emission at 50% of the peak (see Table 1 in Paper I). 
The value of -Nch^cn is ~ 5 times lower than the col- 
umn density derived by means of the RD method (see 
Sect. I4~2~T1) 
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As already mentioned for G31.41+0.31, this 



Fig. 24. Serendipitous detections toward the center of 
G31.41+0.31 obtained with the PdBL The vertical labels 
along the profiles point out the identified lines, reported 
in Tables ^2 and El and the unidentified peaks, listed 
in Table El The y-axis scale is expressed in brightness 
temperature. The dotted vertical lines are corresponding 
numbers indicate the position of the CH3CN (12-11) K- 
components in the upper part part of each spectra, and 
of the 13 CH3CN (12-11) if-components in the lower part 
part. 



is due to the fact that the models do not take into account 
the dumpiness of the region. 

The dynamical mass computed with these values as- 
suming that the toroid is seen edge-on by the observer is 
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Mdyn ~ 34 M . This value is much smaller than the mass 
of the core, and, as for G31.41+0.31, suggests that such a 
structure is unstable against gravitational collapse. The 
infall of material onto the embedded millimeter source 
is also suggested by the line widths, AV", measured in 
the toroid (see Fig. EB). Note that for G24.78+0.08 the 
line widths have been measured from CH3CN, because 
it has not been possible to properly measure them from 
CH 3 13 CN, probably due to an opacity effect. As seen in 
the figure, AV steadily increases toward the center of the 
toroid and peaks at the position of the millimeter contin- 
uum source. The estimate of v- m { derived from the differ- 
ence between AV at the position of the 1.4 mm continuum 
source and that at the edge of the toroid (see Sect. 13. 2. "ty 
is Vjjai ~ 1.2 km s -1 in this case. This value is slightly 
lower than the u ro t obtained from the models or from the 
velocity gradient in the toroid. The accretion rate derived 
for such a V[ n f is M acc ~9x 1CP 3 M Q yr _1 . Such a huge 
accretion rate supports the disk- accretion scenario also for 
the formation of the massive YSO embedded in core Al 
in G24.78+0.08. 

5. Other molecules in G31.41+0.31 and 
G24.78+0.08 

Figures [21 and |2U show some portions of the spectra ob- 
tained toward G31. 41+0.31 and core A2 in G24.78+0.08, 
respectively. Taking into account that the spectra for both 
core Al and A2 in G24.78+0.08 are quite similar, we 
have only plotted the spectra toward A2, where the de- 
tections are more clearly visible. In these figures one can 
see that besides the expected CH3CN and 13 CH 3 CN pat- 
terns, numerous other lines have been serendipitously de- 
tected. These lines have been identified, when possible, 
by using the Jet Propulsion Laboratory, Cologne, and 
Lovas databases for molecular spectroscopy. Emission due 
to C2H5CN, C 2 34 S, HNCO, DC0 2 H, and CH3OCHO-A 
high excitation (100-270 K) transitions have been mea- 
sured. For these lines, Table ^2 lists the transition, the 
rest frequency, the upper level excitation E u (K), and line 
parameters obtained with a Gaussian fit, i.e. the velocity 
Vlsr (km s" 1 ), the FWHM line width (km s" 1 ), the peak 
brightness temperature Tb (K), and the integral under the 
line jT B dV (K km s" 1 ). 

A special case is represented by a peak at 220.46 
GHz (see Table 1110 . where the only reasonable candi- 
date to our knowledge is the transition Jk_k + =20 2 .is- 
19 3 ,i 7 of glycolaldehyde (CH 2 (OH)CHO), whose excita- 
tion (120 K) well fits with the temperature measured to- 
ward G31.41+0.31 and G24.78+0.08. Glycolaldehyde is an 
isomer to methyl formate (CH3OCHO), here detected at 
220.45 GHz. Whereas the detection of CH3OCHO is not 
surprising since it is a common species in hot cores (e.g. 
van Dishoeck & Blake ESP, CH 2 (OH)CHO has been de- 
tected up to date only toward Sagittarius B2N (Hollis et 
al. I2000[) . The reasons could be that complex interstellar 
molecules seem to prefer the C-O-C backbones with re- 
spect to the C-C-0 ones (see Fig. 2 of Hollis et aL lMTUjl . 



Table 13. Unidentified lines detected toward 
G31.41+0.31 and G24.78+0.08 



Frequency a 
(MHz) 


Spectral resolution 
(MHz) 


Source 


220665.91 ±0.15 


2.5 


G31, G24 


220564.07 ± 0.29 


1.3 


G31 


220551.97 ±0.41 


1.3 


G31 


220426.40 ±0.18 


2.5 


G31 


220361.37 ±0.50 


2.5 


G31, G24 


220308.68 ± 0.43 


2.5 


G31 


220258.91 ±0.31 


2.5 


G31 



(a) The frequencies have been measured from the spectra 
toward G31.41+0.31, due to a better signal-to-noise ratio of 
the profiles, assuming a systemic velocity Vlsr= 97 km s _1 . 



Note that other CH2(OH)CHO lines lie inside the spec- 
tral windows here investigated, but unfortunately they are 
closely blended with other bright lines. In conclusion, we 
can only claim a tentative detection of CH2(OH)CHO to- 
ward G31.41+0.31 and G24.78±0.08, which calls for fur- 
ther observations to confirm the occurrence of glycolalde- 
hyde emission in hot cores. Note that there is a difference 
between the line rest frequency and the frequency of the 
observed peak (~2 MHz), as indicated by a slightly differ- 
ent peak velocity from the Vlsr of the region ( Table ITT|) . 
If the detection of glycolaldehyde is confirmed, this would 
imply that the rest frequency listed in the spectral cata- 
logues should be refined. 

Table El reports the observed frequency, the corre- 
sponding spectral resolution and the line candidates with 
their frequency and excitation for those spectral features 
that could not be unambiguously identified. The candi- 
date transitions have frequencies that fall inside the spec- 
tral features, and were selected on the basis of the fol- 
lowing criteria: (i) emission lines known to be hot core 
tracers, such as CH3OCH O-E, C 2 H 3 CN, an d C 2 H 5 OH 
(e.g. van Dishoeck & BlakeESEl Cazaux et al. 2003), and 
(ii) lines with excitations in agreement with the temper- 
atures measured toward G31.41+0.31 and G24.78±0.08 
(E u <; 500 K). Finally, Table Il3l reports the peak frequen- 
cies and the relative spectral resolution of the unidentified 
lines toward G31.41+0.31 and G24.78±0.08. In conclu- 
sion, the gas associated with G31.41±0.31 results more 
chemically rich than that in G24.78±0.08: this could be 
caused by a different evolutionary stage or it could simply 
be a luminosity effect being G31.41±0.31 more luminous 
than G24.78±0.08. 

In order to understand which regions are traced by 
these complex molecules, the maps of the emission due to 
the unblended lines is shown in Figs. [21)1 (for G31.41±0.31) 
andE3(for G24.78±0.08). The distributions of C 2 H 5 CN, 
HNCO, CH3OCHO-A, and CH 2 (OH)CHO (greyscale) 
have been compared with that of the K=0, 1, and 2 
components of CH 3 CN (contours). For G31.41+0.31, it 
is possible to see that the emission due to these molecules 
comes from the central region where the protostars are lo- 
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Table 11. Serendipitous detections toward G31.41+0.31 and G24.78+0.08 (Al and A2) 



Molecule 


Transition 


Rest frequency 


E u 


Vlsr 


FWHM 


T B 


J Tb dV 






(MHz) 


(K) 


(km s" 1 ) 


(km s _1 ) 


(K) 


(K km s" 1 ) 


G31.41+0.31 


C2H5CN 


(25 2i24 -24 2i23 ) 


220660.92 


143 


96.3 ±0.1 


7.7 ± 0.2 


50.6 ± 3.3 


415 ± 12 


34 nc 


(24 23 -23 2 3j 


220645.47 


187 


96.1 ± 0.5 


7.2 ± 0.5 


1 T A _l_ A O 

17.4 ± 4.2 


133 ± 15 


HINOO 


(10i,9-9i,s) 


220585. 20 e 


102 


9o.o ± (J. 4 


12.2 ± 1.2 


on T J_ A 

o9.7 it 4.3 


518 ± 39 




fin. , q~ „ \ 

(LU(,,i »6,3 J 


220561.16 s 


139 




c A -L 1 c; 
0.4 it 1.0 


14 44-9 1 

l^.^t I Z.l 




CH 2 (OH)CHO h 


(20 2 ,i 8 -19 3 ,i7) 


220463.87 


120 


93.1 ±0.8 h 


8.2 ± 1.8 


22.1 ± 3.7 


193±36 


CH3OCHO-A 


(25ii,i 5 -26 9 ,i 8 ) 


220445.79 


272 


96.0 ± 1.1 


7.9 ±2.4 


10.9 ±3.4 


92±26 






G24.78+0.08 


core Al 








C 2 H 5 CN b 


(25 2 , 24 -24 2 , 23 ) 


220660.92 


143 


111.4 ±0.1 


4.3 ±0.2 


24.3 ± 2.0 


111 ± 5 


HNCO d 


(10i, 9 -9i, 8 ) 


220585. 20 e 


102 


111.8 ± 1.7 


7.0 ± 1.7° 


19.9 ±1.7 


148 ± 28 


DC0 2 H f 


(10 6 ,4-9 6 ,3) 


220561.16 s 


139 


111.7 ±0.2 


3.9 ±0.7 


8.9 ± 1.1 


37±5 


CH 2 (OH)CHO h 


(20 2 ,1 8 -19 3 ,17) 


220463.87 


120 






3.0 ±0.5' 


41±3 ! 






G24.78+0.08 


core A 2 








C 2 H 5 CN b 


(25 2 , 24 -24 2 , 23 ) 


220660.92 


143 


110.4 ±0.1 


6.7 ±0.4 


33.4 ± 1.4 


237 ± 11 


c 2 34 s c 


(24 23 -23 23 ) 


220645.47 


187 


110.5 ±0.6 


6.1 ± 1.0 


10.7 ± 1.6 


70 ± 13 


HNCO d 


(10i, 9 -9i, 8 ) 


220585. 20 e 


102 


110.9 ±0.1 


7.5 ± 1.2 C 


35.3 ±0.8 


281 ± 5 


DC0 2 H f 


(10 6 ,4-9 6 ,3) 


220561.16 s 


139 


109.7 ±0.2 


7.2 ±0.4 


16.1 ± 1.1 


124±7 


CH 2 (OH)CHO h 


(20 2 ,i8-19 3 ,i7) 


220463.87 


120 


107.1 ±0.3' 


6.4 ±0.7 


3.9 ±0.4 


26±3 



(a) Spectra taken toward the peak positions of G31.41±0.31 and cores Al and A2 in G24.78±0.08. 

(b) Blended with the 220665.91 unidentified line (see Table HH. 

(c) Blended with the K =5 component of CH3CN (12-11). 

(d) Blended with the K=6 and K=4 components of CH3CN and 13 CH 3 CN (12-11), respectively. 

(e) The HNCO (10i, 9 -9i,8) line splits into six hyperfine components with AF=0,±1 spread over 0.6 MHz and thus not resolved 
with the present spectral resolution. The frequency refers to the F-F =10-9 line. 

(f) This line could be blended with the high excitation (E n = 625 K) C 2 H 3 CN (10 3 ,7-10 2 , 8 ) line at 220651.39 MHz. 

(g) The rest frequency corresponds also to the DC0 2 H (106,5-96,4) line, which has the same excitation. 

(h) According to the databases for molecular spectroscopy, glycolaldehyde results to be the most probable candidate to identify 
the present emission at this frequency. If confirmed, this would be the first detection of CH 2 (OH)CHO outside the Galactic 
center (see text). 

(i) Bad fit because of non-Gaussian profile. 



Table 12. Other molecular tracers detected toward G31.41±0.31 and G24.78±0.08 



Obs. frequency a 


dv 


Source 


Line candidates 




(MHz) 


(MHz) 




Transition 


Frequency (MHz) 


E U (K) 


Blended lines 


220702.17 ± 0.15 b 


0.7 


G31 


CH3OCHO-E (104,7-93,7) 


220702.75 


43 








CH3OCHO-E (24 3 , 22 -24 2 , 23 ) 


220701.52 


169 








CH3OCHO-E (24 3 , 22 -24i, 24 ) 


220701.36 


169 








CH 2 (OH)CHO (40 9 ,32-40 8 ,33) 


220701.03 


511 


220498.34 ± 0.08 c 


1.3 


G31 


C 2 H 5 OH (19 5 ,i4-19 4 ,i5) 


220500.70 


250 








CHj 3 CH 2 CN (74,3-63,4) 


220498.40 


29 


220345.51 ±0.16 


2.5 


G31, G24 


C 2 H 5 OH (7 3 , 4 -7i, 7 ) 


220346.86 


91 








C 2 H 3 CN (9 4 , 5 -10 3 , 8 ) 


220344.45 


55 



(a) The observed frequencies have been measured from the spectra toward G31.41±0.31 due to a better signal-to-noise ratio of 
the profiles. 

(b) Blended with a non Gaussian profile, which produces a plateau in the observed profiles (see Fig. 1241 . A possible candidate 
for the emission is C 2 H 3 CN b (10 3 ,7-10 2 , 8 ) at a frequency of 220695.42 MHz. 

(c) This line could be blended with the very high excitation (E u = 1005 K) CH3OCHO-E (46o,46 -45 4 ,4i) line at 220667.14 
MHz. 



cated as indicated by the positions of the 1.4 mm (white 
cross) and 7 mm (white dots) peaks. This result con- 
firms the presence of a hot molecular core associated 
with G31.41±0.31, in agreement with the increase of tem- 
perature and density toward the center of the region as 



measured with CH3CN (see Fig. [TJ). Moreover, Fig. I2"u1 
shows that the emission from different tracers peaks 
at slightly different positions, indicating that C2H5CN, 
HNCO, CH3OCHO-A, and CH 2 (OH)CHO are tracing dif- 
ferent portions of the hot core. Figure [Z7\ shows a simi- 
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G31 .41 +0.31 




a (J2000) a (J2000) 

Fig. 26. Overlay of the PdBI CH3CN (12-11) emission averaged under the K = 0,1, and 2 components (contours) and 
the C2H5CN (25 2 ,24-24 2 , 23 ) emission (top left), the HNCO (10i )9 -9i, 8 ) emission (top right), the CH 2 (OH)CHO (20 2 ,i 8 - 
193,17) emission (bottom left), and the CH3OCHO-A (2511,15-269.18) emission (bottom right) toward G31.41+0.31. 
Contours are the same as in Fig. EJ Greyscale levels range from 0.10 to 1.18 Jy beam -1 in steps of 0.18 Jy beam -1 for 
C 2 H 5 CN, from 0.10 to 0.94 Jy beam -1 in steps of 0.14 Jy beam -1 for HNCO, from 0.10 to 0.46 Jy beam -1 in steps of 
0.12 Jy beam -1 for CH 2 (OH)CHO, and from 0.01 to 0.13 Jy beam -1 in steps of 0.02 Jy beam -1 for CH3OCHO-A. The 
synthesized beam is shown in the lower left-hand corner. The white cross marks the position of the 1.4 mm continuum 
emission peak, and the white dots the 7 mm continuum emission peaks detected by Hofner (private communication). 



lar scenario for G24.78+0.08: the C 2 H 5 CN, HNCO, and 
CH 2 (OH)CHO (CH3OCHO-A has not been detected to- 
ward G24.78+0.08) emission traces the two hot cores, Al 
and A2, confirming that each of them is associated with 
a distinct site of star formation and is hosting YSOs. 

6. Conclusions 

We have analyzed the millimeter high angular resolution 
data, obtained with the BIMA and IRAM PdBI interfer- 
ometers, of the dust and gas emission toward the high- 
mass star forming regions G31.41+0.31 and G24. 78+0.08. 
The aim was a thorough study of the structure and kine- 
matics of the rotating toroids detected in both hot cores 
in Paper I. 

The CH 3 CN (12-11) emission of the toroids in 
G31. 41+0.31 and core Al in G24.78+0.08 has been mod- 
eled assuming that it arises from a disk-like structure 
seen edge-on, with a radial velocity field. For G31.41+0.31 
the model properly fits the data for a velocity v mt ~ 



1.7 km s 1 at the outer radius i? ut — 13400 AU and 
an inner radius i?i nn — 1340 AU, while for core Al in 
G24.78+0.08 the best fit is obtained for u rot ~ 2.0 km s -1 
at R out ~ 7700 AU and R inn ~ 2300 AU. Unlike 
the Keplerian disks detected around less luminous stars 
such as IRAS 20126+4104 ( e.g. C esaroni et al. H9"97|l . 
NGC 7538S (Sandell et al. I27TO . or M17 (Chini et 
al. l2004|l . these toroids are not undergoing Keplerian rota- 
tion. From the modeling itself, however, it is not possible 
to distinguish between constant rotation or constant angu- 
lar velocity, since both velocity fields suitably fit the data. 
The best fit models have been computed assuming for the 
toroids a temperature gradient of the type T oc i? -3 / 4 , 
typical of accretion circumstellar disks (Natta 2000 ) , with 
a temperature at the outer radius T out ~ 100 K for both 
cores. 

The Mdyn needed for equilibrium derived from the 
models is much smaller than the mass of the cores. This 
suggests, as already pointed out in Paper I, that such 
toroids are unstable and undergoing gravitational col- 
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C 2 H 5 CN 





HNCO 
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12 s 3 
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Fig. 27. Overlay of the PdBI CH3CN (12-11) emission averaged under the K = 0,1, and 2 components (contours) and 
the C 2 H 5 CN (25 2 ,24-24 2 ,23) emission (top left), the HNCO (10i )9 -9i, 8 ) emission (top right), and the CH 2 (OH)CHO 
(202,i8"193,i7) emission (bottom left) toward cores Al and A2 in G24.78+0.08. Contours levels range from 0.10 
to 1.00 Jy beam -1 in steps of 0.18 Jy beam -1 . Greyscale levels range from 0.10 to 0.70 Jy beam -1 in steps of 
0.10 Jy beam -1 for C 2 H 5 CN, from 0.10 to 0.60 Jy beam -1 in steps of 0.10 Jy beam -1 for HNCO, and from 0.02 to 
0.10 Jy beam -1 in steps of 0.02 Jy beam -1 for CH2(OH)CHO. The synthesized beam is shown in the lower left-hand 
corner. The white crosses denote the position of the 1.4 mm continuum emission peaks. The triangle marks the position 
of the VLA 1.3 cm continuum emission peak (Codella et al. [1997 ) . 



lapse. The collapse is supported also by the fact that the 
CH 3 13 CN or CH 3 CN line width measured in the cores in- 
creases toward the center of the toroids in G31.41+0.31 
and G24.78+0.08 Al. The estimates of v^f and M acc are 

- 2 km s -1 and - 3x 10 -2 M yr -1 for G31.41+0.31, and 

- 1.2 km s -1 and - 9 x 10 -3 M yr -1 for G24.78+0.08 
Al. Accretion rates that large could weaken the effect of 
stellar winds and radiation pressure thus allowing further 
accretion on the star. This would support theories accord- 
ing to which high-mass stars form like their lower mass 
counterparts through non-spherical accretion with large 
accretion rates. 

The values of T ro t and -/Vch 3 cn, derived by means 
of the RD method, for both G31.41+0.31 and the sum 
of cores Al and A2 (core A of Codella et al. 1997) 
in G24.78+0.08 are in the range 132-164 K and 2- 
8 x 10 16 cm -2 . 

For G31.41+0.31, the most plausible explanation for 
the apparent toroidal morphology seen in the lower K 
transitions of CH3CN (12-11) in Paper I is self-absorption, 
which is caused by the high optical depth and the existence 



of a temperature gradient in the core. Indeed, the -/Vch 3 cn 
presents a peak toward the apparent dip in the CH3CN 
brightness temperature map. 
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Fig. 25. Serendipitous detections toward the peak posi- 
tion of core A2 in G24.78+0.08 obtained with the PdBI. 
The vertical labels along the profiles point out the identi- 
fied lines, reported in Tables ITT1 and ITS! and the uniden- 
tified peaks, listed in Table The y-axis scale is ex- 
pressed in brightness temperature. The dotted vertical 
lines are corresponding numbers indicate the position of 
the CH3CN (12-11) if-components in the upper part 
part of each spectra, and of the 13 CH 3 CN (12-11) K- 
components in the lower part part. 



